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1. Introduction 



Following 't Hooft's pioneering work [l], instanton physics has played an im- 
portant role in understanding the dynamics of asymptotically-free theories. Mostly, 
explicit instanton calculations were limited to the semi-classical approximation. A 
two-loop calculation which is particularly relevant for supersymmetry (SUSY) can be 
found in ref. Q. 

There is no consistent regularization method which preserve SUSY ||. In this 
sense, SUSY is similar to a chiral symmetry. Whether SUSY is a true symmetry at 
the quantum level, or not, can be decided only by non-perturbative studies. This 
question was previously addressed in the semi-classical approximation, where it was 
found that SUSY is preserved. Important works are based on the assumption that 
this remains true beyond the semi-classical approximation, and that SUSY is an exact 
symmetry in the continuum limit [4-10]. 

We begin with the investigation of Ward identities in the most general path 
integral setting. A Ward identity reads 

{SO) = (SS O) + (Sfi O) + ( Hilbert-space surface terms ) , (1.1) 

where the bold letter S denotes any infinitesimal transformation. Besides the familiar 



terms containing the variations of the action (SS) and the measure (Sfi), eq. fll.l]) 
contains also (the expectation value of) Hilbert-space surface terms, which arise from 
integration by parts in the functional integral. 

In more detail, in each non-perturbative sector the path integration is over a 
given set of collective coordinates ( n that characterize the classical background, as 
well as over the amplitudes of the quantum-field modes. The expectation value of a 
(gauge invariant) operator O is 

(O) = J d n ( e - s ^- Wl - W2 (0) c . (1.2) 

Here S c \ is the classical action, and the semi-classical measure, exp(— W±), is the 
product of functional determinants and a tree-level jacobian [lj]. W2 is the sum of all 
connected bubble diagrams, and (O)^ is the sum of all diagrams with external legs 
defined by the operator O. All diagrams are obtained here by integrating over the 
quantum fields, keeping the collective coordinates fixed. Eq. ( |TT2] ) applies to both 
bare and renormalized quantities. 

Eq. (|i~T| ) is derived in Sec. 2. Restricting our attention to the SUSY variation, 
which we denote 5, we write SS and SfM as spectral sums, and show that they both 
vanish mode by mode. Thus, (SO) may fail to be zero only due to the Hilbert-space 
surface term. Explicitly 

(SO) = fd n (^- (e-**-^-** (S( n 0> f ) , (1.3) 

where S( n is the SUSY variation of the n-th collective coordinate. 
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The tools of functional differentiation and integration used in Sec. 2 are very 
useful in clarifying the algebraic structure that underlies eq. ( |1.3| ). However, manipu- 
lations of the (formal) path integral should be supported by a detailed diagrammatic 
calculation. In Sec. 3 we turn to the instanton sector of four-dimensional SUSY gauge 
theories. Only the integral of the total p-derivative is not automatically zero (p is 
the instanton size). The anomalous term may thus be represented as an amplitude 
involving a zero-size instanton. An anomaly due to a topological singularity was pre- 
viously found in supersymmetric quantum mechanics in the presence of fluxons with 



a fractional magnetic flux |11 



We consider the SUSY Ward identity (|1.3| ) at the leading non-trivial order for a 
specific bi-local operator in super Yang-Mills theory. Our main result is the anomalous 
Ward identity eq. ( [3.10|) , and its local version eq. ( |3.38|) . Both equations are valid 
in any renormalization prescription where perturbation theory is supersymmetric. 
While the leading contribution to the p — > surface term on the r.h.s. of eq. (|3.10|) is 
a tree diagram, a straightforward calculation of the expectation value of 50 involves 
loop diagrams. The proof of eq. ( |3.10| ) thus requires careful separation of short- and 
long-distance contributions. We also generalize the result to supersymmetric QCD. 
An important question is how eq. ( |3.10| ) is modified by higher-order corrections. A 
renormalization-group argument indicates that the anomalous term should remain 
finite (and non-zero) after the inclusion of next-order logarithmic corrections. 

The limit of a vanishing instanton size, p —>■ 0, is studied in more detail in Sec. 4. 
In correlation functions involving operators of sufficiently high dimension, the lower 
limit of the p-integral may diverge. Under these circumstances one should perform 
non-perturbative subtractions. Our main result is that no non-perturbative subtrac- 
tion can possibly modify eq. ( |3.10|) . Therefore eq. (|3.10|) constitutes a supersymmetry 
anomaly. We also discuss how the anomaly should arise on the lattice. 

There are important non-perturbative sectors where one is not expanding around 
an exact solution of the classical field equations. This is true in the one-instanton 
sector of theories with a scalar (Higgs) VEV, as well as in the instanton-antiinstanton 
sector. Aiming to cover these cases too, we have developed Landau-gauge Feynman 
rules for any non-perturbative sector of a renormalizable gauge theory (Appendix B). 
The Feynman rules define a systematic expansion provided the considered correlation 
function is infra-red convergent, and provided the background field is "almost" a 
classical solution. 

Affirming the validity of the perturbative expansion must be done on a case 
by case basis. (A more detailed discussion of the one-instanton sector of SUSY- 
Higgs theories will be given elsewhere.) However, the basic Ward identity ( |1.3|) can 
be derived in a completely general way using the Feynman rules of Appendix B. 
This is done in Appendix C. We give a fully detailed diagrammatic derivation in the 
(physically less interesting) case of a theory without a gauge symmetry. Because of 
its technical complexity, we outline the generalization to gauge theories but omit the 
arithmetical details. 

Much of the technical material contained in this paper is devoted to various 
detailed proofs. A first acquaintance with the algebraic structure may be obtained by 
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reading Sec. 2 up to subsection 2.2. In order to understand the content of eq. ([3.10D it 
is enough to read Sec. 3 up to subsection 3.1. The local form of eq. ( |3.10| ) is discussed 
in subsection 3.5, and the renormalization-group argument is given in subsection 3.6. 
The proof that eq. ( |3.1(J| ) cannot be modified by subtractions is given in subsection 4.1. 

The notation used throughout the paper is defined in Appendix A. In Appendix D 
we explain how the general formalism of Sec. 2 works in the case of translation 
invariance. Our conclusions are given in Sec. 5. 



2. Path-integral derivation of SUSY Ward identities 

In this section we derive the anomalous Ward identity ( |1.3|) using the tools of 
functional integration and differentiation. While the (unregularized) path integral 
is a formal construct, the derivation serves to clarify both the underlying algebraic 
structure, and the role of physical boundary conditions. The results of this section 
are supported by the detailed diagrammatic calculations of Sec. 3 and Appendix C. 

The notation of Appendix A is used below. The general and the SUSY Ward 
identities (eqs. ( |1.1| ) and ( |1.3j ) respectively) will be derived using the transverse-field 
path integral defined in the first four subsections of Appendix B. (The rest of Ap- 
pendix B, which is devoted to the construction of the tree- level propagators etc., is 
not needed for this section.) 

2.1. The field variation 

The independent variables of the path integral include the collective coordinates 
( n , the gauge degrees of freedom u a (x), and the amplitudes of the quantum modes 
f3 p and ^ p where 

fc*) = t X B p {x)h , = £ 4{x)% , (2.1) 

•tp •'p 

and Xp( x ) (Xp( x )) is the eigenfunction corresponding to /3 P (\& p ). In this paper the 
quantum bosonic field, $(x), is assumed to be transverse. Namely, it obeys the back- 
ground gauge condition Vft a f3{x) = (eq. ( P-5Q ) in addition to the orthogonality con- 
ditions (& ;m |/3) = (eq. ( |B.16p ) where b- m is the covariant derivative of the classical 
field with respect to ( m . The above constraints may be summarized as {b. M \(3) = 0, 
where the capital index m runs over the ordinary collective coordinates and over the 
infinitely-many collective coordinates associated with local gauge transformations, see 
Appendix B.4. 

We begin by expanding an arbitrary variation of an elementary field in terms of 
variations of the independent path-integral variables. For the variation of a boson 
field, SB(x), one has 

SB(x) = 6$(x) + (b. n (x) + p. n (x)) 6Cn + (n a (b(x)) - igT a p(x)) 6cu a (x) . (2.2) 
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Covariant (^-derivatives of the quantum field are obtained by differentiating the eigen- 
functions (see also Appendix B.3) 



= £ x£„(*) & , (2-3) 
while the quantum-field variation is, by definition 

6$(x) = i^^(x)6$ p . (2.4) 

The field variation acts on the quantum amplitudes (3 P and not on the eigenfunc- 
tions. The differential operator Q a generates infinitesimal gauge transformations of 
the classical bosonic field (Appendix B.l). The variation of a fermion field, 8^f(x), is 
expanded as 

6V(x) = dif(x) + 5( n % n (x) - ig 8u a (x) T a , (2.5) 

where 

*in(x) = t Xp-M % > = i Xp (x) 8% , (2.6) 

in analogy with the bosonic case. (The reader should note that 8^(x) and 8^(x) 
have different meanings.) For the SUSY variation, ordering matters in eq. ( |2.5| ). It 
follows from eq. fl2.4|) that S/3(x) obeys the same orthogonality constraints as does 
/3(x). This allows us to extract the variations 8( n and 8u a (x). Considering eq. (|2.2| ) 
and taking the inner product (|A.5| ) with b- n leads to 

5( m = C-l l {(h n \ry)+ i g(blT a $\5u a )} , (2.7) 



where C mn is defined in eq. (B.14) and we have used eq. (P-7| ). We have specialized 
to the SUSY case, where the bosons' variation is linear in the fermions, see eq. ([A. 71 ). 
(For any other variation, simply replace 5 — > 8 and T^/ — > 5B.) Note that b- n is 
0(l/g), and d~( m is 0(g). We next apply on both sides of eq. (|2.2j ). Using the 
gauge-fixing constraint ( |B.5| ) and eliminating 5( m we find 

8u a = (^- 1 ) ab \n t ^ ri) 

+z^(^- 1 ) ab |6| m T fe /3) Cj, (fe ;n |r#) . (2.8) 
Here = JFq" 1 — T^T-^Tq + • ■ •, where JF and .F int occur respectively in the 



tree-level and interaction ghost lagrangians (eq. (]B.37| ) ) . Eqs. (|2.7|) and fl2.8|) are 
summarized by the compact formula (cf. eq. ( |B.31| )) 

5C M = C^ N (b, N \T^) . (2.9) 

The variation of a given quantum amplitude is obtained by taking the inner 



product of eq. (|2.2| ) or fl2.5|) with the corresponding eigenmode. Explicitly, 

5% = (x? |r*) - 8C M (xf |Am) , (2.10) 



for bosons, and 

6* p =(4\6*)-6Cu(4\%u)- (2-11) 

For fermions. The last term is these equations compensates for the (^-dependence of 
the eigenmodes. 

2.2. Ward identities 

The variation of any local operator can be expressed in terms of the variations of 
the independent variables. One has 

50 = TO, (2.12) 
where T is the functional differentiation operator 

r = p^ + pc,i-, (2.i3) 



and q p stands for any quantum-field amplitude, cf. Appendix B.3. In practice we 
will be interested in gauge invariant operators, which effectively restricts the last 
sum to n = n only. The expectation value of SO is given by the functional integral 
(eq. ( F33D ) 

r Det C 

(SO) = / d n ( Vq i e~ s TO . (2.14) 
J Det 2 C 

Integrating by parts we obtain the generic Ward identity for gauge invariant operators 



(so) = (ss o) + (Sn o)+ I d n c^-\fv q e ' s s (nO 

J d( n [J Det^Co 



(2.15) 



where SS and 8 fx are respectively the variations of the action and the measure. Ex- 
plicitly 8S = TS and 

_^rio g {^} + £^ + ^. (2,6) 

iDet^CoJ V d 1p 9 Cn 



The last term on the r.h.s. of eq. ( |2.15| ) is the Hilbert-space surface term of eq. ( |1.1|) . 



The corresponding surface term for the quantum amplitudes q p is zero because of the 
Gaussian integration. 

In the rest of this section we consider SUSY Ward identities. We compute SS 
and 8fi, and prove that they both vanish mode by mode. This leaves us with the last 
term in eq. fl2.15|) . When the functional integration over the quantum fields is done, 



this term is recognized as the r.h.s. of eq. ([L 

The SUSY transformation is off-diagonal in that it maps bosons into fermions 
and vice versa. In view of this, the invariance of the measure, Sll = 0, is an expected 
result. Nevertheless, the below derivation gives us valuable information on how the 
cancelation of the various contributions to Sll works in practice. We also find that 
the physical boundary conditions play a direct role. The same elements reappear in 
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the diagrammatic proof of Appendix C which is more rigorous and, at the same time, 
technically more involved. 

The presence of a Hilbert-space surface term in Ward identities is not an unfa- 
miliar situation. In appendix D we apply the generic Ward identity ( j2.15 ) to the 



case of translation invariance. We show that momentum is conversed because the 
Hilbert-space surface term coincides in this case with a spacetime surface term. 

2.3. Variation of the action 

The functional variation of the action consists of a volume integral and a surface 
integral. The latter may in general depend on the choice of functional variables, and 
it is our objective to determine it. We start with the action principle 

J d d x (SC - d^(5Q n„)) = J d d x 5Q (-n M;M + dC/dQ) . (2.17) 



As usual C = C(Q,Q tfJj ) and ILj = dC/dQ^. Substituting eqs. (|2.2|) and (|2.5|) into 
the r.h.s. of eq. ( |2.17|) and integrating by parts leads to 

J d d x 5Q (-n^ + OC/dQ) = J d d x 5q (-U^ + OC/dQ) (2. 18a) 

+5Q n J d d x (Q^IL, + Q, n {dC/dQ)) (2.18b) 
-SCnfda^j;. (2.18c) 

Gauge invariance of the action was used. (Note that dC/dQ = dC/dq.) Here 

j; = Q ;n n M . (2.i9) 

Let us write T = T q + 7^ in correspondence with the two terms on the r.h.s. of 
eq. (|2.13| ). Since Cn-derivatives act on the classical field or on the wave-functions 
(eq. ( |B.19| )), expression ( |2.18b| ) is equal to T^S. Now, in terms of the mode ampli- 
tudes, the bilinear part of the action reads 

* (2) = ~|>P<g, (2.20) 

where X p is an eigenvalue of the small fluctuations operator, Lx P = \Xp- Conse- 
quently 

T q SW = t\ p 5q p q p =(5q\L\q). (2.21) 

In infinite volume L must act on q(x), and not on 5q(x). The expression on the r.h.s. 
of row (|2.18a|) is therefore equal to T q S. Putting everything together we arrive at the 



following result 

T3 = J d d x 8C- i da^ (5Q n M ) + 6( n j j; . (2.22) 
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Eq. ( |2.22|) is completely general. In the SUSY case, 5C is a total derivative and 

TS=i da, S, + 8( n I da, J% , (2.23) 



where S, is the SUSY current. In order to understand the last surface term we employ 
a finite volume cutoff. This term is then completely determined by the boundary 
conditions and, for any boundary conditions which ensure the hermiticity of the small 
fluctuation operators, 

^tZ(7 M Xp,«n M = 0. (2.24) 
We note that eq. ( 2. 24] ) holds for each eigenmode Xp separately. This implies the 



operator statement § da, J™ = 0. 

2.4. Variation of the measure 

Computing the variation of the measure is a straightforward algebraic task. The 
first term on the r.h.s. of eq. ( 2.16|) is 



Tl0 A^A =ClKTC^-kCo)llT(CoU- (2.25) 
The classical-valued matrix Co depends only on the collective coordinates, and 

T {Co)kl = ^Cm {Co)kL;M — 3(n (Co)jf L;n ■ (2.26) 

The last equality is true since (Cq) KL;M = for any collective coordinate ( M which 
corresponds to an exact symmetry. The latter always include (local) gauge transfor- 
mations and (global) translations. Hence the non-zero terms always correspond to a 
subset of the ordinary collective coordinates ( n . 

In the calculation of the field-dependent part TC\ we use the bosonic closure 
relation 

^t\x?)(x?\+±\K)(b M \, (2.27) 

~ —1/2 

where b M = (C ) MN b- N are normalized. Including the contribution from C one has 

- (b. KiL \r*) - 5( M {h K \B, M , L ) . (2.28) 

Next, using eq. Q2.9| ) we have 

= -e^c KL;n 8c L + c^(b ;K \rif) n . (2.29) 

What remains is the contribution from the functional trace. For the fermions 

t ^Sr 1 = ^ t - c - (M r **) > ( 2 - 3 °) 



where eq. ( 2.11 ) and the closure relation for fermions were used. For the bosons, 



starting from eq. ( |2.10| ) one has 

t ^ = -*C t {x B p \x B p ;n) - KM t { X B p \h) {h K ; M \ X B p ) (2.31) 

Using the bosonic closure relation ( p.27|) we find after a few algebraic steps 



£^ = -5Cn£(x P 1x^)-^£(L|&Mj-^nC-iC KL;n 
Ofjp Jp^ 

+\Kn {C^) LK {C^) KL , n + KuC'l (b, K \B. L;U ) . (2.32) 
Collecting all terms we obtain 

^ = $(n £(-) F (Xp|Xp;n) + Kn £ (b M \b M - n ) 

+6( M Cl\ (b, K \B, M , L - B. L;M ) . (2.33) 

The first term on the r.h.s. is a spectral trace over (transverse) bosons and fermions. 
As with the variation of the action we now employ a finite volume cutoff. The modes 
are normalized, namely, (xp|Xp) — 1- Hence (x P |Xp;n) = 0. Similarly, (b M \b M ;n) = 
0. Finally, as we show below, the commutator (the last term in eq. ( |2.33| )) is in fact 
a spectral trace over the ghosts field, which is zero for the same reason. Thus, <5/i 
vanishes mode by mode. 

2.5. Ghosts contribution 

We last consider the commutator in eq. ( [2.331 ). Using Appendices B.3 and B.4 
and inserting a complete set of ghosts eigenstates, the commutator can be rewritten 

as 

(b. K \B. M . L -B. L;M ) = (b. K \{tt a -igT a p)F* L )=i^C Kp F pML , (2.34) 
where F^ L is a generalized field strength and 

F pML =(c a p \F a ML ). (2.35) 

Hence 

C IA - (b;K | -B;M;L -B;L;A/) = V ] C ' LK C Kp F pML = V ] Fp Mp , (2.36) 

where explicitly 

F pMp — ~ ( c pI c p;a/) • (2.37) 



We see that the commutator in eq. (|2.33j ) is a spectral trace over the ghost field, 
which vanishes in finite volume too. 

In summary, using the (formal) tools of functional differentiation and integra- 
tion we showed that the only anomalous term in SUSY Ward identities is a surface 
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term in Hilbert space, cf. eq. (|1.3|) . The contributions of SS and Sp (cf. eq. (|2.15|) ) 
vanish mode by mode. We note that the classical equation of motion was nowhere 
used. Consequently, eq. ( |1.3| ) should hold in any non-perturbative sector, regardless 
of whether or not one is expanding around an exact classical solution. A detailed 
diagrammatic proof of this statement is given in Appendix C. 



3. Anomalous supersymmetry Ward identities 

In this section we turn to the one-instanton sector of super Yang-Mills (SYM) 
theory. Our main result is a SUSY Ward identity whose anomalous term is non-zero 
and unambiguous. The same anomalous term is obtained using any renormalization 
scheme where perturbation theory is supersymmetric. Moreover, as shown in Sec. 4, 
no non-perturbative subtraction can modify this result. Similar statements apply to 
supersymmetric QCD (SQCD). 

We begin (subsection 3.1) with a leading-order calculation of (SO) for a specific 
gauge-invariant bi- local operator in SYM (eq. ( |3.10| )). Eq. ( |1.3|) is used in the calcula- 
tion. As already mentioned in the introduction, the result arises from (the p-integral 
of) a total p-derivative. 

In the next two subsections we confirm eq. ( |3.1(J| ) by a detailed diagrammatic 
calculation. Working at fixed p we show (subsection 3.2) that every diagram which 
contributes to (50) p is related to a diagram which contributes to the p-derivative of 
(8pO) p . While (5pO) p is a tree diagram to leading order, one encounters loop dia- 
grams in a direct same-order calculation of (SO). In subsection 3.3 we prove that 
the (tree-diagram) r.h.s. is equal to the (one- loop renormalized) l.h.s. of eq. ( |3.10| ). 



A key role is played by the recursion relations of Appendix B.6. (While this dis- 
cussion is given in terms of a concrete example, it actually proves the validity after 
renormalization of eq. (|1.3j ), whenever its r.h.s. consists of tree diagrams only.) 

In subsection 3.4 we generalize the result to SQCD, and in subsection 3.5 we 
discuss the local form of the anomaly. In subsection 3.6 we give a renormalization- 
group argument that the anomalous term should remain finite, and non-zero, after 
the inclusion of higher-order logarithmic corrections. 

We now give a quick review of (super) Yang-Mills instantons taking the gauge 
group to be SU(2). In SYM one expands around an exact, scale-invariant classical 
solution, and the Feynman rules are considerably simpler than those of Appendix B. 
The (regular gauge) instanton field is 



c _ 2 r} CfJjU (x xp) u ^ 
g [x — xq) z + p A 



where f\ clLV is the 't Hooft symbol [jl|. (fj CfMU is antisymmetric in the last two indices; 
Vcfiu = Ccnv if p,v = 1,2,3, and fj c4:U = 5 CU .) There are eight gauge-field zero modes 
(b- n in our generic notation). Their explicit (unnormalized) form is a c . As explained 
in Appendix B.l, each zero mode constitutes of an ordinary derivative of the classical 
field with respect to a collective coordinate, plus a compensating (proper) gauge 
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transformation to enforce the background gauge a M;n = 0. We now list the various 
zero modes. First, for the dilatation zero mode, covariant and ordinary differentiations 
coincide i.e. af„„ = af, „. For the zero modes associated with the translation collective 
coordinates Xq, the compensating gauge transformation is generated by uj c v = a c v . 
Explicitly a^ v = a^ v — a v = f^. 

Last we consider the three isospin zero modes a„. a /. In a singular gauge, one has 

o^ = £^a;J, (3.2) 

where „ , „ 

_ % (x - x Q f 5 a a , p 2 \ 

The last term in parenthesis generates a proper gauge transformation, whereas the 
first term generates a global SU(2) transformation. During the quantization proce- 
dure described in Appendix B.2 one introduces global SU(2) collective coordinates 
( a > alongside with the translation and dilatation ones. As follows from eq. (|3.3|) , 
the isospin modes obey a^ a i = a M)0 / — D^uj a ' where gu%, = 5%, p 2 /{(x — xq) 2 + p 2 ), 



in agreement with the general formula eq. ( B.6 ). Then, in gauge invariant correla- 
tion functions the isospin collective coordinates factor out, yielding a group- volume 
factor 0. 

We next turn to the gaugino. Because of the chiral nature of instanton amplitudes 
it is often more natural to use a Weyl notation, where the Majorana- field gaugino is 
split as A ma j orana — > X a , A A . Similarly, we will distinguish between the SUSY variations 
5 a and 5 & (both of which are accounted for in eq. ( |A.7| )). In SU(2) SYM there are 
two pairs of gaugino zero modes. In normalized form they are given by 

* SS (PT a = — T( %-^2- /? = 1,2, (3.4) 

7T ({X - Xq) 2 + p 2 Y 

A {a)a ~ n ((x-x y + p iy ' a " 1;2 - (3 ' 5) 

In terms of Weyl fields, the instanton-sector fermion propagator [ITS] obeys 




o ) v ; v p p( x >y) ) 

(3.6) 

where Pp is the projector on the (left-handed) fermionic zero modes, and (•••)□ refers 
to the tree-level instanton propagators. In the more compact Majorana notation the 
same equation reads 

L F G F (x,y) = 5\x,y)-P F (x,y). (3.7) 
The transverse vector-field propagator |12| obeys 



L B G B (x, y) = 5\x, y) - P B (x, y) - P« (x, y) , (3.8) 



where P B projects on the eight bosonic zero modes (cf. eq. QB.58 )) and P" is the 



longitudinal projector ( |B7BTj ). (In SYM Q -> £> M ; also, P B L B = P" L B = 0, cf. 

eq. mm .) 
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Last we give the semi-classical measure (d n ( exp(— S c \ — W\) of eq. (|1.2| )), which 
reads 

2 10 rr 6 g- 8 A 6 1 dpp 3 d A x , (3.9) 

where Ai is the one-loop renormalization-group invariant scale. In SYM the functional 
determinants, cf. eq. (|B.45| ), cancel each other exactly. The matrix Co (eq. ( |B.14a|) ) is 



diagonal. The entries related to the translation and the dilatation modes have mass 
dimension zero, while those related to the three isospin zero modes have dimension 
minus two. The p 3 factor in eq. ( |3.9|) thus arises from the isospin entries of det^Co- 
In (potentially) singular cases, the integration over the translation collective co- 
ordinates Xq should always be done before the p-integration. This guarantees transla- 



tion invariance. Also, it implies that the integral of the total xo-derivative in eq. ( L3j) 
vanishes. The isospin collective coordinates too do not give rise to any Hilbert-space 
surface term, because group- integration is compact. 

3.1. Super Yang-Mills 



In eq. ( |1.3| ), a non-zero result may arise in the instanton sector (only) from the 
p-integral of the total p-derivative. In this subsection we compute the anomalous 
Ward identity 



6 a Of (0) XX(y)) = Vi^Wrlhu (3 1Q) 

Here 5 a acts on everything to its right, AA = A^ X aa and 

Of = (D; b (e bcd X ca ■■<% X*)) ((&yf\tP) , (3.11) 

where 

D, = d,- igT c A c lx = Dp — igT c a^ , (3.12) 
and a 1 = \ f\i tlv <J iiV . Using eq. (|1.3|) , the l.h.s. of eq. (|3.10|) is expressed as a limit 

(6 a Of(0)XX(y)) = 2Vg~ 8 AfJdp-?-{ P 3 ((5 a p)Of(0)XX(y)\ 

= -2Vg- 8 Af lim p 3 ((6 a p) Of (0) AA(y)) . (3.13) 

For p 2 ^> y 2 the p-integrand behaves like y/p 7 (see subsection 4.1 below). The integral 
is infra-red convergent, and the p = oo boundary point drops out from eq. ( p,13| ). 



Eq. ( 3.13 ) says that the anomalous term can be expressed solely in terms of an 
amplitude of a zero-size instanton. This mathematical statement has one misleading 
aspect. Namely, the dominant contribution to the p-integral which leads to eq. ( |3.10|) 



comes from finite-size instantons with p 2 ~ y 2 . Eq. ( ^.10[ ) is insensitive to the contri- 
butions of vanishingly- small instantons. If we restrict the integration to p < p < oo 
for some < p 2 y 2 , the relative change in the result will be 0(p 2 /y 2 ). The limit 
p — » is smooth. We return to these observations in Sec. 4. 



12 



With eq. ( |3.13|) at hand the calculation is straightforward. We begin with 5 a p 



(cf. eq. (|2.7p). To leading order one has 



lQn 

where 5 a A* = (i/V2)M a6t A cA , and 

16tt 2 



<W = I d ' x ( d </ d P) > ( 3 - 14 ) 



g 2 



d 4 x(da c Jdp) 2 , (3.15) 



is the pp-entry of C . Note that 5 a p = 0(g) since da^/dp = 0(1/ g). 

As can be seen from the above equations, 5 Q p contains a X a (x) field. To arrive at 
eq. (|3.10| ), one gaugino field from XX(y) is contracted with X a (x) to form a propagator 
Gp(y,x) (cf. eq. ( |3.7| )), and the second one is saturated by a X sc (a) mode. (If A a 
is contracted with one of the A's in Ojf(0), one obtains a contribution that vanishes 
in the limit p — > due to an extra damping factor of p 2 /y 2 -) Now, because of the 
operator OJf, the integration over the instanton position is dominated by xq ~ p. 
The integral in eq. Q3.14 ), in turn, is dominated by x ~ y. In the calculation of 



the p — > boundary term, this allows us to replace the (singular gauge) fermion 
propagator Gp(y,x) by a free propagator. We find 

(»MM >lOTi = |^i^ , (3.16) 

where we have indicated which (superconformal) Grassmann amplitude has been 
integrated over. 

The other three zero modes saturate the operator OJf. The xo-integration yields 



(At this order, only the background covariant derivative contributes, cf. eq. Q3.ll 



We now put together eqs. (|3.16| ) and (|3.17| ), include a symmetry factor e Q/3 , and 



substitute everything on the r.h.s. of eq. ( |3.13j ). The result is eq. ( |3.10| ) 
3.2. Diagrams 



In this subsection we describe in detail how eq. (|3.10| ) works at the level of Feyn 



man diagrams. This provides a concrete leading-order example of the general dia- 
grammatic identities of Appendix C. The renormalization of the loop diagrams en- 
countered in the derivation will be discussed in the next subsection. 

Taking a slightly different course from Appendix C, our starting point is the 
observation that, in SYM, 

r d A xdJs^(x)0) =0, (3.18) 
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for any (multi) local operator O. Here, O will be taken to be the bi-local operator in 
eq. ( ^.10[ ) . In the generic notation of Appendix A, SUSY invariance of the action can 
be expressed as the off-shell identity 

d ll S fl = (p ] L B + S B )Ti>+m(L F i> + S F ), (3.19) 

where 

5f = ^, 5b = ^ + ^, (3.20) 

dm dp dp 

and the interaction lagrangians are defined in eqs. (|B.43|) and ( |B.44| ). The expressions 
inside the parenthesis in eq. (|3.19| ) are recognized as the (bosonic and fermionic) equa- 



tions of motion. Below we will repeatedly use eqs. (|3.18| ) and (|3.19|) . (We have used 
the notation d^S^ in the above off-shell relation to distinguish it from the operator 
d^S^ the latter is discussed in subsection 3.5 below.) 

The first diagrams that contribute to (5 a O) are 0(1/ g 2 ) in comparison to the 
r.h.s. of eq. ( |3.10| ). At this order, (S a O) vanishes for rather trivial reasons, as we now 
explain. For any SUSY theory, the leading-order terms in eq. ( |3.19|) are 

d^S® = Sgr® + M (0) L F V , (3.21) 



where eqs. ( |3.20| ) and (P-44j ) were used. We have expanded 

5V = 5¥ 0) + <^ (1) + c^ (2) , (3.22) 

and the superscript counts how many quantum fields occur in each term. In SYM, 
eq. (|3.21|) reduces to 



d,S^ = 6 a \® p p . p ¥ = i/(2V2) 0%/,., p . A* , (3.23) 



Like eq. (|3.19|) , eq. (|3.23| ) is an off-shell relation in which only the classical equation 
of motion is used. The fermion propagator emanating from d^S^ ends on any one of 
the five A's in the operator O. The other four A's are saturated by zero modes. We 



next eliminate this fermionic propagator using eq. ( |3.7|) . Only the delta-function from 
eq. ( |3.7| ) contributes to (S a O). The terms with the projector P F cancel each other 
by antisymmetry. These steps are depicted in Fig. 2. (See Fig. 1 for our Feynman 
rules.) To this order, the diagrams contributing to (S a O) involve the fermionic zero 
modes and the classical field, and nothing else. As in eq. ( |3.18|) , the vanishing result 



is obtained before the integration over any of the collective coordinates. 

We now turn to the more interesting next order, that corresponding to eq. ( 3.10| ). 



First, there is a set of disconnected diagrams which is trivially zero, since it consists 
of the same diagrams considered above (Fig. 2) times a sum of bubble diagrams. 
(The bubble diagrams (Fig. 3) were discussed in detail in ref. 0. Note that, besides 
the familiar two-loop diagrams, Fig. 3 contains a one-loop diagram coming from the 
off-diagonal terms of the jacobian (|B.13|) or, equivalently, from the interaction term 
on the last row of eq. ( |B.39| ). This is the only place where that term occurs at this 
order.) 
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We next consider Fig. 4. The shown diagrams are related by interchanging the 
insertion of d^S^) with one of the zero modes. Thanks to this antisymmetrization, 



the fermionic projector terms of eq. (|3.7|) cancel out. After the application of eq. ( |3.7| ). 
the second diagram gives a contribution to SO. 

In the first diagram of Fig. 4, the propagator emanating from d^S^ is attached 
to a fermionic interaction vertex. On the l.h.s. of Fig. 5 we show only the connected 
part of this diagram containing d^Sj®'. Now, the terms in eq. ( |3.19| ) which are linear 
in both the fermion and the boson fields read 

d^S® = ft L B ™ + M (1) L F ^ + W (0) S F . (3.24) 

After the application of eq. ( |3.7| ), the l.h.s. of Fig. 5 gives rise to an insertion of (minus 
the integral of) the last term in eq. (|3.24j) (middle Fig. 5). Using eq. (|3.18| ), the latter 
is traded with an insertion of the first two terms on the r.h.s. of eq. ( [3.241 ). Explicitly, 
this insertion is 

(r\^|L B |/3) + . (3.25) 

This step is depicted on the r.h.s. of Fig. 5. (Eq. ( |3.25|) coincides with eq. ( |C.2| ) up 
to the replacement 5^ — > 5^^.) 

We now apply eqs. (|3.8|) and (|3.7| ) to the diagram containing insertion (|3.25| ). 



The terms with the fermionic projector cancel by antisymmetry as before. All other 
terms are shown in Fig. 6. The last diagram, which involves a longitudinal projector, 
will be discussed later. The first diagram on the r.h.s. of Fig. 6 is a contribution to 
50. (This contribution is a tree or a one-loop diagram, depending on whether the 
two external legs correspond to different spacetime points or to the same one.) The 
second diagram vanishes by the same Fiertz rearrangement used in proving the SUSY 
invariance of the action. (Here we are ignoring the need for regularization, see the 
next subsection). 

The third diagram on the r.h.s. of Fig. 6 is a contribution to the r.h.s. of eq. ( |3.10| ). 



The bosonic zero mode b- n at the upper-left corner of the triangle, the A" at the same 
spacetime point and the thick dashed line representing Cq , together form the leading 
order of 5 a ( n . The other bosonic zero mode is attached to a fermion line emanating 
from a fermionic zero mode. Denoting the zero mode by x i this part of the diagram 
gives the (^-derivative = —Gp (L F ). n x F0 - 

We now consider two more diagrams with insertion ( |3.25| ) . The third diagram on 



the r.h.s. of Fig. 7 gives the (^-derivative of the fermionic propagator emanating from 
5( n . Fig. 8 contains a one-loop tadpole which is actually zero in SYM. In general, the 
third diagram on the r.h.s. of Fig. 8 gives the (^-derivative of the (logarithm of the) 
functional determinants. 

Let us consider separately the boson, fermion and ghost contributions in the 
second diagram on the r.h.s. of Fig. 8. The fermion-loop contribution is by itself zero 
after a Fiertz rearrangement. The ghost-loop diagram will be considered later. For 
the sum of the remaining one-(bosonic)-loop diagram, plus the first diagram on the 
r.h.s. of Fig. 7, we use 

dpS® = S$g (b) /W> A + W^L F * + W {l) S F . (3.26) 
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This allows us to trade the two bosonic-loop diagrams with a single diagram contain- 
ing an insertion of S^^Lp Using eq. ( |3.7| ) once more, we obtain the contribution(s) 
to 50 comine; from 5\W (Fig. 9). 

Another type of diagrams with insertion ( 3.25|) is obtained by attaching the 



bosonic propagator to a quantum gauge field coming from a covariant derivative 
(cf. eq. (3.12[)). One obtains two contributions to 50 and a contribution to the r.h.s. 



of eq. (ft.lOp , in which the classical field in the background covariant derivative is 
differentiated with respect to ( n (Fig. 10). 

We next attach the bosonic propagator from insertion ( |3.25| ) to a vertex coming 



from the expansion of the jacobian (Fig. 11). The second term on the r.h.s. of Fig. 11 
gives the (^-derivative of the bosonic zero mode contained in 5( n itself. 

The third diagram on the r.h.s. of Fig. 11 gives the ^-derivative of logdet^Q). 
Particularly relevant for eq. ( |3.10| ) are p-derivatives. As explained earlier, the p- 
dependence of det Cq comes from the entries related to the three isospin zero modes. 
In the matrix C\ (eq. ( |B.14b| )), the isospin-isospin entries (Ci) a > c > involve the integral 
of ge a bctic' a 1,a' ( c f- eo ±- After taking the product with da® /dp, the result can 

be written as ig tr a Mi(X / [a) c > , da^/dp}. This yields the p-derivative of det 5 Co, since 



d d 

— — CLu c' — TT" 

op op 



<V, C ' = ^- Dp u) d = -ig 



dp 



+ D^. (3.27) 



The contribution of the last term is zero after integration by parts (which is allowed 
since duj c i/dp vanishes rapidly enough at infinity) using D^a^ a i = 0. (The above is 
an example of the commutator formulae of appendix B.3.) 

Last we discuss all diagrams with the longitudinal projector. Compare the last 
two diagrams in Fig. 6. These diagrams have the same topology, but the meaning of 
the various elements is different. Recalling the spectral decomposition of the ghost 
propagator, the analog of the zero mode b ]n is now Q a Cp where is a ghost eigenstate 
(see below eq. ( |3.8| ), and Appendix B.4). Similarly, the corresponding inverse eigen- 
value (associated with the ghost propagator) is the analog of Cq 1 (the thick dashed 
line). Thus, those elements which constitute 5( n in the third diagram, correspond in 
the last diagram to 5uj p = (cp\5uj a ) (both to leading order). 

Similarly to what we did with the bosonic modes b. n , the insertion of Vt a c p on 
a fermion line can be traded with a local gauge transformation with parameter c^, 
acting on the field(s) at the end(s) of the line. The sum of the resulting diagrams, 
where the gauge transformation acts on the fields of O, vanish by gauge invariance 
of O. In addition there are diagrams that vanish thanks to gauge invariance of the 
semi-classical jacobian, or of 5uj p itself. One last piece is provided by the ghost- 
loop diagram that has remained from the second term on the r.h.s. of Fig. 8. It is 
recognized as a contribution to the gauge transformation of the product c a p F^f 
occurring in 8ui p . 

The above completes the diagrammatic analysis of eq. ( |3.10|) , exception countert- 
erm diagrams. These will be discussed in the next subsection. 
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3.3. One-loop renormalization 



Several types of loop diagrams occur in the calculation of eq. ( |3.10| ), and the 



corresponding divergences are renormalized by counterterms. In this subsection we 
complete the proof of eq. (|3.10|) , assuming that renormalized perturbation theory is 
super symmetric. (This statement means that the perturbative S-matrix is super- 
symmetric, and that perturbative matrix elements of composite operators fall into 
supermultiplets.) We assume that the counterterms were constructed using the back- 



ground field method . For definiteness we will refer to dimensional regularization, 
but the discussion generalizes to any other consistent regularization as well. As men- 
tioned earlier, the below arguments are actually sufficient to prove eq. ( |1.3| ) after 
renormalization, whenever its r.h.s. consists of tree diagrams only. 

As explained in the previous subsections, the two-loop bubble diagrams that occur 
on the l.h.s. of eq. ( |3.10| ) are multiplied by zero. Thus, we need not concern ourselves 
here with the corresponding counterterm diagrams. (There are delicate points in the 
renormalization of the bubble diagrams, see ref. [0]; these will become relevant at 
higher orders.) 

We next consider the ghost-loop diagram contained in the second term on the 
r.h.s. of Fig. 8, and all the (loop) diagrams with the longitudinal projector P". As 
discussed in the previous subsection, the sum of these diagrams is zero by gauge 
invariance. Since (background) gauge invariance is preserved by the regularization, 
this cancelation continues to hold. 

The remaining divergences arise from one-loop diagrams containing boson or 
fermion propagators, and no ghost propagator. Renormalization of these diagrams 
requires us to face two obstacles. First, dimensional regularization treats a loop as a 



single entity. Equations like eq. (|B.64|) , ( p.8|) or (|3.7|) do not hold for d — 4 ^ 0. We 



must therefore replace these equations by ones that hold for arbitrary d. The second 
(related) complication is that SUSY, as expressed by eq. ( |3.19| ), is broken by terms 
proportional to d — 4. 

We now explain how dimensional regularization is applied to an instanton-sector 
diagram. First, repeating n times the recursion relation for the generic fermion prop- 
agator G F , eq. ( |B.68|) , we obtain 

G F = (l - P£ xact ) G^ c J2 {-V F G^ c ) k + G F (-V F G^ c ) n+1 . (3.28) 

fc=0 

(The superscript "vac" denotes free- field quantities, and V F — L F — Vf c . Notice the 
terms with the projector pj, xact ; which "cut open" any loop containing G F .) We now 
take n (finite and) high enough, such that any loop containing L F times the last term 
in eq. (|3.28|) will be finite. Eq. (P.64j ) may be applied to the last term, and the result 



is 

n-l 



L F G F = (l™ c + V F ) G™ c Y, (-VfGf") + L^ C GT (-V F G™) n - P£ xact . (3.29) 



k=0 



We have used (see below eq. ( |B.68| )) 



pexact (_y F Q^y = pexact _ (3 3^ 
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Eqs. ( |3.28|) and ( 3.29| ) provide the key for regularization. Comparing eq. ( ft .29 ) 
with eq. ( |B.64| ), we see that the term with pj, xact has been separated out explicitly. 
The remaining terms, which replace the delta-function in eq. (|B.64j ), are identical to 
what one would have found in a standard perturbative expansion. Any (loop) diagram 
containing (only) these perturbative-expansion terms is amenable to dimensional (or 
any other) regularization. (If L™ c G™ c (x, y) = 5 4 (x — y) is true in the regularized 
theory, then eq. ( |3.29| ) reduces to eq. (p.64| ).) 

(In SYM, the bosonic and fermionic propagators in the instanton sector obey 
eqs. ( |3.8| ) and ( |3.7| ) respectively, and the above analysis is applicable to both of them. 
In the bosonic equivalent of eq. (|3.29|) , pj, xact is replaced by Pb + P "■ The ghost 
propagator contained in the longitudinal projector (cf. eq. ( |B.61| )) obeys the standard 
recursion relation ( [B.66|) . In the most general case Gb is defined by eq. (|B.54|) . The 
last term in this equation cuts open any loop since (like projectors) it involves a 
localized source. As for Gb, it is convenient (cf. Appendix B.6) to consider first 
£ 7^ 0. Then Gb obeys eq. (|B.48| ) and the standard recursion relation (|B.66|) . Using 
eq. ( |B.59| ) , all term involving l/£ may be replaced by expressions that have a smooth 
£ -> limit.) 

Let us now examine, in comparison, the background-field fermion propagator in 
the vacuum sector. This propagator obeys the standard recursion relation (|B.66| ) 
and, therefore, identities analogous to eqs. ( ft.28| ) and ( |3.29| ) but without projector 
terms. Similarly, the vacuum-sector identities for the background-field transverse 
boson propagator contain the longitudinal projector, but no analog of the projector 
Pb- 

When the manipulations based on eq. ( ft,19|) are carried out in the vacuum sector, 
all (finite) violations of SUSY that survive the removal of the regularization can be 



canceled by an appropriate (non-supersymmetric) set of counterterms p4|| . The can- 
celation of SUSY violations arising from the loop regularization continues to hold in 
the instanton sector. As a result, whenever LbGb or LpGp arise inside some (ampu- 
tated) one- loop diagram (from the application of eq. ( |3.19| )) this yields a renormalized 
diagrammatic identity in the instanton sector which differs from the corresponding 
one in the vacuum sector precisely by the projector terms pj, xact and Pb ■ As shown in 
the previous subsection, all diagrams with an insertion of pj, xact cancel by antisym- 
metry, while those containing Pb combine with diagrams that involve the functional 
jacobian to form the total ( n - (and in particular p-) derivatives. 

We now explain how the last statement works in practice. Consider first the 
loop diagram on the l.h.s. of Fig. 7, which originates from the self-energy diagram 
in Fig. 12(a). There is a corresponding counterterm diagram, shown in Fig. 12(b). 
(As mentioned earlier, the sum of one- loop tadpoles in Fig. 8 is zero in SYM.) Now, 
the r.h.s. of Figs. 7 and 8 also contains loop diagrams, to which we have to apply 
eq. ( ft.26| ), cf. Fig. 9. That equation, which is an expression of SUSY of the action, 
is broken in dimensional regularization by an amount proportional to d — 4. When 
multiplied by id — 4) _1 coming from the loop divergence, a (finite) explicit breaking 
of SUSY may result. This explicit breaking is canceled, however, by an appropriate 
non-supersymmetric set of counterterms. Consequently, the renormalized equalities 
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represented by figs. 7, 8 and 9 hold after adding the counterterm diagram. 

A second set of one-loop diagrams can be found in Figs. 4, 5 and 6 (provided 
both of the external legs go to the same spacetime point). The divergences in these 
figures correspond to the composite operators A A or OJf (eq. ( p.llj )), or to their SUSY 
variations. For example, 5 (XX) = (i/y/2) cr^F^X. However, in general 

5[AA]«^(z/V2)[a M ^A]«, (3.31) 

where [• ■ •] W is the corresponding one-loop counterterm. The l.h.s. of inequality ( [3.31D 
is shown in Fig. 13(a), while its r.h.s. is Fig. 13(b). (The counterterms are normally 
0(g 2 ), but when they involve a classical field they become 0(g).) Again, the countert- 
erms are designed to compensate for any discrepancy that may have arisen from the 
loop regularization. With the counterterm diagrams added, the equalities of Figs. 4, 
5 and 6 hold after renormalization. Finally, the loops of Fig. 10 are treated in a sim- 
ilar manner. We leave it for the reader to work out the corresponding counterterm 
diagrams. This completes the proof of eq. ( |3.1U| ). 

At higher orders, the renormalization of eq. (|3.10j ) is more complicated, because 



one must deal with the divergences arising from the coupling between the discrete- 
and the continuous- index parts of the jacobian. This problem was addressed in ref. @. 
A related problem is that one must deal with the divergences of 5 p. In subsection 3.6 
below we use a renormalization-group argument to determine the next-order logarith- 
mic corrections to eq. ( |3.10|) . 



3.4. SUSY theories with matter 

It is easy to generalize eq. ( |3.10| ) to SU(2) SQCD. Let the number of flavors be 



Nf. (By convention, for SU(2) each flavor corresponds to two chiral supermultiplets 
in the fundamental representation. Each (massless) Weyl-fermion field has one zero 
mode.) One has 



N t \ „19 5„ i»6-iV 



S a Of (0) XX(y) n i&(0)) = ~ 2 " y ^t|^k^ , (3.32) 
M>k(q) = [ d*xe i<lx ^ h (x) . (3.33) 



where 



We have assumed y 2 m 2 <C 1 where m is a generic matter-field mass. This implies that 
only the zero modes contribute in eq. fl3.33|) . As shown in subsection 4.1 below, both 
eq. ( |3.10|) and eq. ( |3.32|) cannot be modified by subtractions and, hence, constitute a 



supersymmetry anomaly. Eq. (|3.32| ) reflects a general feature, namely, the anomalous 
term may be a polynomial in some of the external momenta (here the polynomial in 
g M reduces to a constant). 

For SU(N), One should replace Ojf by an operator whose generic structure 
is X(D fl X)(D 2 X)(XX) N ~ 2 . The generalization of eq. (|3.10|) will have on its r.h.s. 
Cjv g l ~ m A^ Nf times similar spacetime dependent factors. We have not worked 
out the numerical constant Cjv- 
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3.5. Local form of the anomaly 



The local version of eq. ( |1.3| ) is obtained by applying the SUSY variation of the 
fields only at the point x = z. To this end, one multiplies the l.h.s. of eqs. (|2.2| ) 
and ( j2.5| ) by 5 4 (x — z). One can now repeat the construction of Sec. 2. An easy 
way to read the local variations from those of Sec. 2.1 is to associate the S 4 (x — z) 
factor with the matrix V. Whenever Y^> occurs inside some integral, it is replaced 
by the corresponding integrand at x = z. For the local variations of the collective 
coordinates one obtains 

5( m (z) = C-l \b\ n {z) Tif(z) +igj d A xb] n (x) T a p(x) 5u a (x; *)} , (3.34) 



where 



5u a (x;z) = T- b \x,z)Q b \z)T^(z) 

+ig f d 4 y^- b \x,y)bUy)T b P(y) C£bUz)T$(z). (3.35) 



The variation of the measure is zero as before, and the local form of eq. ( |1.3| ) reads 
(S^z) O) = contact terms - / ' d n ( J- ( e -Sa-W-W a (^(^ )J ? ( 3 . 36 ) 

where the contact terms are the expected variations of the (multi)local operator O. 
The last term in eq. ( |3.36| ) is, by definition, a matrix element of the operator d p S p {z). 
Another way of reaching these results is to keep track of the diagrammatic identities 
of Sec. 3.2, but now without performing the integral in eq. ( |3.18| ). For the local form 
of eq. (|3.10| ) we find 

d v (S ua (z) Of (0) XX(y)) = contact terms - (d v S va (z) Of (0) \\{y)) , (3.37) 
where the matrix element of d v S ua (z) is 

dAa(z) Of (0) AA(y)) = -2Vg- 8 Al lim p 3 ((8 a p(z)) Of (0) XX(y) 

2 u - 5 Af r] hpv z v (y - z) K rjbxa y a _ , _ _ ]f3 . . 

~ 7gt (z^y {{y-zfY (y2)2 V t ,Acr p cx K a x a T \ a . [S.SX) 

Eq. ( |3.38| ) is consistent with locality of d u S ua (z). As expected, one recovers the r.h.s. 
of eq. ( |3.10| ) by integrating over z in eq. ( |3.38| ). 

3.6. Renormalization group considerations 

At the next-to-leading order, the Hilbert-space surface term on the last row of 
eq. ( |3.13| ) may pick up logarithmic corrections. Depending on the power of log p, the 
result of the p — > limit could blow up, remain finite, or vanish. 

In fact, the logarithmic corrections may be either logp or logy 2 . This is because 
the limit in eq. (|3.13|) is not uniform: while the effective range of the xo-integration 
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scales to zero with p, the distance scale y 2 is kept fixed. Taking this delicate point 
into consideration, we now give a renormalization-group (RG) argument that the log p 
factors arising from 5 a p and Off should cancel those arising from two-loop bubble 
diagrams. 

Given a renormalization point p, we rewrite eq. ( |3.10|) to one higher order as 



g\p) (6* Of (0) XX(y)) = -^g\p) 2 ^ 2 ^f^ (1 + logs) . (3.39) 

The two-loop RG-invariant scale is 

= p 6 g~\p) exp(-87r 2 /^ 2 (/i)) . (3.40) 



The reason for including the g 5 (p) factor in eq. ( |3.39|) is as follows. First, the op 



erator 5 a p is RG-invariant because the collective coordinate p is independent of the 
renormalization point p, and the SUSY variation 5 a respects RG-invariance. More 
generally, the basic building blocks of RG-invariant operators are gF^ v and gX (or 
gX). We have therefore multiplied Off (0) and XX(y) by appropriate powers of g(p). 

We now analyze the expected log p corrections, starting with the contribution of 
the two- loop bubble diagrams . The semi-classical instanton measure ( |3.9| ) involves 
the product A^g'^^p). RG invariance requires that, for size-p instantons, the loga- 
rithms arising from the bubble diagrams should turn this product into h\g~ A (p~ l ). 
Next consider the operator S a p. As can be seen from eqs. fl3.14|) and (|3.15|) , the cou- 
pling constant explicitly contained in S a p is associated with the classical field. The 
logarithmic corrections to 5 a p should modify that g(p) to g(p~ l ). 

Finally, RG invariance of g 3 Ojf means that the exponential involving its anoma- 
lous dimension, exp[/jfe ' dg' jo(g') I (3(g')], must be proportional to g 3 (p^ 1 )/ g 3 (p). 
In other words, the renormalization of Off provides log p factors which are just enough 
to turn g 3 (p) into g 3 (p~ 1 ). Putting together the expected logp factors arising from 
all sources we find that they cancel each other. (Of course, it is important to confirm 
this conclusion by a direct next-order calculation.) Moreover, all higher-order correc- 
tions to both jo{g) an d (3(g) cannot give rise to logp terms in the solution of the RG 
equation. 

Finally, if the scale y 2 is varied, one expects the renormalization of XX(y) to 
generate log?/ 2 factors that turn g 2 (p) into g 2 (\y\~ l ). We thus expect the next-order 
result to be 

g\p) (5 a Of (0) XX(y)) = -A^M" 1 ) ^^^f^y, {3M) 

Assume now that y e A® 1. While SYM is confining, the contributions of 
sectors with additional instanton-antiinstanton pairs should be damped compared 
to eq. ( ft.lOD (or eq. (|3.41|) ) by extra powers of y 6 Af. What we are calculating 
is thus the leading-order result in an expansion in the physical parameter y®A\ = 
ex V (-8n 2 /g 2 (y^)). 
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4. The p — ► limit 



The limit of a vanishing- size instanton, p — > 0, may be singular in correlation 
functions that involve operators of sufficiently high dimension. In this section we 
investigate several aspects of this limit. Our main result (subsection 4.1) is that the 
anomalous Ward identity eq. ( p,10|) cannot be modified by non-perturbative subtrac- 



tions and, hence, constitutes a supersymmetry anomaly. In subsection 4.2 we discuss 
point splitting, and explain why it does not regularize the operator Off occurring in 
eq. (|3.10p . In subsection 4.3 we describe an ambiguity that arises in a case where 



point splitting can be used and comment on its implications. In subsection 4.4 we 
explain how the anomaly should arise on the lattice. 

4.1. Non-perturbative subtractions 

In this subsection we first examine separately the two correlators ((S a Off) XX) 
and (Off (S a XX)), whose difference is eq. ( |3.10| ). We show that the corresponding 
p-integrals are convergent in the limit p — > 0. 

For comparison, we next consider another Ward identity where the corresponding 
integrals are divergent for p — > 0. In that case non-perturbative subtractions are 
needed, and may in fact be used to recover SUSY. Having worked out this explicit 
example, we list the general properties of non-perturbative subtractions. It then easily 
follows that eq. (ft.lO|) cannot be modified by any non-perturbative subtraction. 



In subsection 3.1 we have computed the correlation function ((5 a p) Off(0) XX(y)) 
(eq. ( p. 13 )) in the limit p — > 0. Let us now generalize the calculation to any p. The 



operator product on the l.h.s. of eq. ( |3.16| ) may be written as 2X sc (y; a) X md (y; a) 
where 

n 2 rr^ a r)n a 

WArV) = | i ^|f. (4.D 
In a regular gauge, the solution is 

bet 

x ind,fe/ \ 99 a l3 , , n \ 

13 {ma) = ¥v{y- Xo y + P 2 - (4 ' 2) 



Performing the Xo-integration we find (ignoring irrelevant constants) 



(4.3) 



where h(s) can be represented as a Feynman-parameters integral (see Appendix E of 
ref. [§) and 

if \ f l + h 1 s + h 2 s 2 + ---, s<l, ( 



Eq. ( p. 10 ) thus involves the integral 



(y 2 f J dp 
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dp-h(p 2 /y"). (4.5) 



This integral is convergent, with the main contribution to it arising from p 2 ~ y 2 . 
For small p 2 /y 2 , the integrand in ( J4.5| ) behaves like py/{y 2 ) 4 - 

We now turn to the correlators ((5 a Ojf(0)) XX(y)) and (0^(0) (5 a XX(y))). We 
claim that both the small- and the large-p behavior of the corresponding integrals is 
the same as in eq. ( |4.5| ) above. In the limit p — > oo this can be established simply on di- 
mensional grounds. Turning to the p — > limit we first examine ((5 a Ojf(0)) XX(y)). 
At tree level, if the product of zero modes at the point y is X sc (y) X (y) then the 
presence of y/{y 2 ) A is evident and, on dimensional grounds, implies a py/(y 2 ) 4 be- 
havior of the p-integrand. Alternatively, if one places both superconformal modes 
at the point y then the factorized ^-independent piece, / d 4 x (5 a Ojf(Q)) xop , is zero 
because the Xo-integral is odd (compare eq. ( |3.17[ )). The first non-zero contribution 
again goes like py/iy 2 ) 4 - At the one-loop level one reaches the same conclusion by 
examining the (small-p behavior of the) fermion propagator in the relevant partial 
wave. Having established this (large- and) small-p behavior of ((8 a OJf(0)) XX(y)), 
eqs. flOD to (|OD imply the same behavior for the p-integrand in (C^(0) (5 a XX(y))). 
(Reaching this conclusion directly is more difficult: the correlator involves a diagram 
with a vector-boson propagator connecting the points and y, and in an instanton 



background this propagator decreases very slowly E2|; see, however, ref. [15 



Before turning to the physical implications of eq. ( |3.10|) , we wish to explain in 



what way things could be different. To this end, we consider another anomalous 
SUSY Ward identity, which is 

2 22 - 5 7r 4 A 6 

(5 & XXF* jk (0) aFX a (y)) = -e ijk gll ^ y2 ^ Ma& V? , ( 4 -6) 



where o~FX a = cr^F^X 013 and 



Ff jk = t a bcF pv F\ p Fl T r} iilv r}j\ p r} k(JT . (4.7) 

The calculation is similar to subsection 3.1, and is in fact simpler. 5&p involves a A 
field, which is saturated by one of the X sc modes, while the other A goes to the 
operator aFX a . Note that eq. (4.6) exhibits a y~ 7 fall-off which is faster than the 



y~ 5 fall-off in eq. fl3.10| ). This kinematic difference will turn out to play an important 
role. 

We now show that the operator XXF^ k requires a non-perturbative subtraction. 
Consider the following correlator 



(XXF? jk (0)XX(y))~e ijk -^[-^, (4 



Al rdp 

sW^ 7 

where on the r.h.s. we have indicated the small-p behavior. The non-perturbative 
logarithmic divergence at small p can be handled as follows. We first restrict the 
integral to p < p < oo where < p. A subtracted operator is defined via 



[XXFf jk ] = XXF? jk - e ijk c 3 log(pp) ^ AA , (4.9) 
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where c 3 is a suitable numerical constant. The subtracted operator [XXFjj k ] yields a 
finite p — > result when substituted into the l.h.s. of eq. (|4.8|). 

Like \\F?j k , the operator 5^ XXFf- k requires a non-perturbative subtraction too. 
We now define the renormalized operator as follows 



[5 & XXFf jk ] = 5 & [XXF* k ] + e ijk 5^ XX . (4.10) 



We have chosen a manifestly non-supersymmetric subtraction. The finite, last term 
in eq. ( |4.10|) was chosen to cancel the r.h.s. of eq. Q4.6|) . Using the vacuum-sector 



result . 

f --F\.(^ rrF\J,A\ = — 

ttV) 

eq. (|4.6|) becomes 



{aFXM crFX a (y)) = -__( ff#l ) o . ^ , (4.11) 



[4 AAF?. fc (0)] o~FX a (y) + lXXF* jk (0)) S & aFX a (y)) = . (4.12) 



Eq. Q4.12p means that SUSY has been recovered in the limit p — > after perform- 



ing the above non-perturbative subtractions. Moreover, by considering additional 
Ward identities, one can verify that the renormalized operators respect the SUSY 
algebra. (Recovering the supermultiplet structure of composite operators by suitable 
subtractions is a reminiscent of the so-called Konishi anomaly [[]].) 

There are two lessons from the above example. First, when non-perturbative 
subtractions play a role, there is at least a chance that SUSY will be recovered. The 
other lesson has to do with the general properties of non-perturbative subtractions. A 
divergence at p — > always arises from a kinematic situation where the instanton sits 
on top of an operator 0{x) of a sufficiently high dimension. Consider the instanton- 
sector correlation function (0(x) Yik ®k{yk)) an d assume that a p — > divergence 
arises only from the operator 0(x). The zero modes and the propagators which are 
functions of yi,y2, • • • , become p-independent in the limit p — > 0. (More generally, 
they are expandable in powers of p; this expansion is used when the leading divergence 
is stronger than log p. ) 

As a result, the dependence on spacetime points of the p — > divergence must be 
that of a vacuum-sector correlation function (0'(x) Y\ k Ok{yk)) va _ c involving another 
operator 0'(x). The operators 0(x) and O f (x) must have the same quantum numbers 
except for their chiral charge, where the mismatch is given by the number of fermionic 
zero modes. Because of the explicit factor which appears in instanton-sector 
correlation functions, the dimension of O' must be smaller than that of O at least 
by 6. (The difference between the dimensions of 0(x) and 0'{x) is 6 (for SU(2)) if 
the divergence is logarithmic; if the divergence goes like some inverse power of p, the 
difference is 6 plus that power.) 

Previously, we showed that there is no small-p divergence in the correlation func- 
tions ((5 a Ojf(0)) XX(y)) and (0^(0) (S a XX(y))). It is now easy to generalize this 
result, and show that no instanton-sector correlation function can have a p — > diver- 
gence associated with the operators OJf or S a OJf '. If such a divergence were to arise 
in some correlation function, then from the spacetime dependence of the divergence 
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one could read off what is the necessary non-perturbative subtraction. However, there 
is no operator that qualifies as a non-perturbative subtraction for Ojf or for 5 a Ojf. 
The mass dimension of Ojf is 7.5. Therefore, for the subtraction one would need an 
operator whose dimension is 7.5 — 6 = 1.5. Evidently, there is no gauge invariant 
operator with this dimension. A similar conclusion applies to 5 a Ojf. Since eq. ( |3.10|) 



cannot be modified by any non-perturbative subtraction, there is no way to recover 
that SUSY Ward identity. Hence, eq. ( |3.10| ) constitutes a super symmetry anomaly. 

One may wonder how the physical implications of eqs. (|3.10|) and (|4.6| ) can be 
so different. Apart from the above mentioned kinematic difference (y~ 7 vs. y~ 5 fall- 
off), there are two more significant differences between the two Ward identities. In 
eq. ( |4.6| ), only the fermionic zero modes and the classical field occur on both sides. As 
shown in ref . [0] , the collective coordinates and the fermionic zero modes constitute a 
finite supersymmetric system. One does not expect that SUSY will be violated when 
these are the only relevant degrees of freedom. On the other hand, in eq. ( |3.10|) the 



l.h.s. involves one-loop diagrams which are outside the scope of the supersymmetric 
calculus of ref. 0, and the potential for an anomaly exists. 

The other qualitative difference between eqs. (|3.10|) and ( |4.6|) is in the corre- 
sponding local Ward identities. In the case of eq. ( |4.6| ), the p — > limit and the 
integration over the point z, where the SUSY variation is performed, do not com- 
mute. Ignoring irrelevant constants, the product X sc \z) da^j 'dp(z) which occurs in 
8^p{z), is p 2 z 2 j[z 2 + p 2 ) 4 . This becomes a delta-function, 5 4 (z), in the limit p — ► 0. 
As a result, the corresponding local matrix element is zero except for z = 0. When 
the only violation of the naive Ward identity is of this type, it is natural to associate 
this effect with a modification of the composite-operator transformation rule, and not 
with a non-zero d^S^(z). 



In contrast, the local form of eq. (|3.10|) (namely Eq. (|3.38|) ) is non-zero for any 
z. Moreover, after the ^-integration one recovers eq. ( |3.10| ), hence the z-integration 
and the p — > limit commute. This result is incompatible with d^S^z) = 0. 

4.2. On point splitting 

In this subsection we address the following question: what happens if one at- 
tempts to define the operators Ojf and o~ a Ojf via point splitting? Since OJf involves 
three gaugino fields, splitting off one of these fields requires the introduction of a par- 
allel transporter to maintain gauge invariance. Our finding are: (a) point splitting 
does not regularize these operators but, rather, introduces new divergences; (b) in 
the instanton sector, the p — ► anomalous term of eq. ( ft.lOD is traded with another 
anomalous term arising from the SUSY variation of the parallel transporter (whose 
evaluation is technically much more complicated). 

For our purpose it is enough to consider the following partial point splitting 

Of[x, e) = (p e »{e bcd X ca (x) o% X^x))) U e f(x, x + e) ((D 2 )^ X^(x + e)) , (4.13) 



where 



U = U(x,x + e) =exp(ig J ds^A^j . (4.14) 
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The SUSY variation of the parallel transporter is 



rx+e / rt \ / rx+t \ 

S a U = igJ dt^expligj ds v A v \ (5 a exp Ug ds v A v \ . (4.15) 



Let us now replace U by 5 a U in eq. Q4.13D . The A field arising from 5 a of eq. ( [4.1 5|) 



may be contracted with any one of the three A fields in eq. ( |4.13| ). For at least one 



of them, the contraction involves same-color fields and, hence, is given in the short- 
distance limit by a free propagator. This free propagator gives rise to a quadratic 
divergence in the t — > x (or t — > x + e) limit at finite e > 0. 

It is well known that point splitting may be used as a (gauge invariant) regular- 
ization in cases such as the chiral anomaly and the Konishi anomaly 0. What is 
common to these examples is that the gauge (and gaugino) fields are external, and 
only matter fields are integrated over. When one integrates over the gauge and the 
gaugino fields too, point splitting ceases to provide a regularization as demonstrated 
above. 

What happens if, nevertheless, one attempts to use the definition ( |4.13| ) in the 
instanton sector with some other method of regularization (e.g. dimensional regu- 
larization)? When the three A's are not at the same point the p-integrand is less 
singular, and behaves qualitatively as 

n 3/(x n \ nW(n A A6 iV^A a ^r}^ yuh( P 2 /y 2 ) f p 2 \ n 

p {(d*p) <V(o, e) x Ky)) p - 7 — \^T7 ) ' ( } 

for some n > 2. The new denominator (e 2 + p 2 )~ n arises from the (possibly differen- 
tiated) wave function of the zero mode at x M = e M . The numerator p 2n corresponds 
to the same wave function in the limit e 2 — ► where eq. ( |4.3|) should be recovered. 
Consequently expression ( Op is replaced by 

Hp 2 /y 2 )(-^) 1 = 0, e 2 >0, (4.17) 




^ + p' 

The p —* surface term thus vanishes for e 2 > 0. In its place, there is now a new 
anomalous term, the one arising from the variation of the parallel transporter ( 4.15| ) . 
After factoring out the leading yj (y 2 ) 3 dependence (and handling the new divergences 
described above!), one is left with an integral over the dimensionless variable p 2 /e 2 . In 
this integral, the parallel tranporters of eq. ( |4.15|) (with A^ taken to be the classical 
field) are 0(1) and may not be approximated by any truncation of their Taylor series. 
The resulting integrals a very complicated and we did not pursue this calculation any 
further. 

4.3. A non-perturbative ambiguity 

If one uses the fermionic equation of motion, the operator Ojf may be traded 
with 

Of = (D; b (e bcd A CQ A#) (e efg . (4.18) 
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After contracting the two epsilon-tensors, one can write Of as a sum of products 

Of = E°IM 2) , (4-19) 

k 

where the (composite) operators O^' and O^ are all gauge invariant. We nay now 
consider the point splitting 



6f(x, x + e) = J2 0^\x) Of\x + e) , (4.20) 

k 



which does not require the introduction of parallel transporters. In the SUSY Ward 
identity, the p — > surface term will vanish again (cf. eq. ( |4.17|) ), and moreover there 



will be no new anomalous terms since there are no parallel transporters. Hence 

, 8 a Of(0,e)XX(y)) = 0, e 2 > . (4.21) 



What is the significance of this result? Let us re-introduce the short-distance 
cutoff p on the p-integral used in subsection 4.1. With both p and e non-zero, one has 

6* Of (0,6) AA(y)) ~ A ' ^ l{fly\ pie) , (4.22) 



where 



W/y\p/t) 



r d 4p{ h{p2/y2) { 



p 2 




(4.23) 



and the approximation p 2 ,e 2 <C y 2 was made on the second row. The result in now 
ambiguous. It depends on the behavior of the ratio p/e in the limit. If we send e — > 
before p — ► 0, eq. ( |3.10| ) is recovered. On the other hand, if we first send p — > and 
only later e — >• 0, then the integrand on the second row behaves like a delta-function, 
S(p), and the result is zero. More generally, if we take the limit with some fixed ratio 
p/e, then p 2 /(e 2 + p 2 ) will interpolate smoothly between and 1. 

We have not been able to resolve this ambiguity in a completely satisfactory way, 
but we have several comments. First, if this problem could not be resolved, this would 
mean an ambiguity in the physical predictions of the theory. The latter conclusion 
seems to us highly unlikely. 

We thus assume that a resolution of the ambiguity should exist. The obvious 
next question is what effect does this have on the SUSY anomaly. Our conclusion is 
that the SUSY anomaly exists regardless of what one does about that ambiguity. We 
will soon argue that, in fact, the prescription of sending e — > after p —>■ violates 
physical principles. Nevertheless, suppose momentarily that that prescription was 
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correct. Eq. ( 4.21 ) would then hold, namely the Ward identity involving % f would 



not be anomalous. However, Ojf and Ojf are related only though the (ultimately 
quantum) equation of motion. The latter is evidently broken by the prescription 
leading to eq. ( f4.21| ). Thus, even if one adopted that prescription, the anomalous 
Ward identity with Off, eq. ( |3.10|) , would still stand, with the conclusion that a 



SUSY anomaly exists. 

In a renormalizable quantum field theory, contributions from the cutoff scale are 
normally associated with some (logarithmic or power-law) divergence. The diver- 
gences are canceled by counterterms, and the ambiguity in the finite parts of the 
counterterms is fixed by renormalization conditions. 

In contrast, the ambiguity we encounter here if of a new type. We find a fi- 
nite contribution of a varying magnitude, coming from an infinitesimal neiborhood 
of the lower end of the p-integral. (In other words, from instantons whose size is 
as small as the short- distance cutoff.) Moreover, these undetermined short-distance 
contributions cannot be attributed to different subtraction schemes, simply because 
(subsection 4.1) they cannot be modified by any subtraction! This state of affairs 
contradicts the principles of renormalization and universality. The only way to avoid 
this problem is to adopt the (unique) prescription where these (otherwise undeter- 
mined) contributions vanish. This prescription amounts to sending e — > before 
p — > (equivalently e/p — > 0). As expected, in this case the Ward identities for Ojf 
and Of agree when the cutoff is removed. 

4.4. SUSY and the lattice 

It is well-known that SUSY is broken by the lattice regularization. Arguments 
that SUSY may be recovered in the continuum limit were put forward in ref. [[16] . 



These arguments are valid in the context of weak-coupling perturbation theory on 
the lattice. However, as we now explain, they do not cover non-perturbative effects. 

In the appropriate momentum (or distance) range, the SUSY anomaly should 
arise on the lattice exactly as in our continuum treatment. The relevant distance 
range is defined by a 2 /y 2 — > where a is the lattice cutoff, and y 6 A^ <C 1 but finite. 
This scaling region, which would be probed in deep-inelastic scattering, is controlled 
by a small coupling constant. The leading dynamical effect in that region is the 
logarithmic evolution of the coupling constant. This is covered by weak-coupling 
perturbation theory (both in the continuum and on the lattice). 

In the scaling region, non-perturbative effects can be studied systematically on the 
lattice just like in the continuum, because they are controlled by the small parameter 
exp(— 87r 2 /(? 2 (|?/|~ 1 )). What needs to be done analytically is to repeat the construction 
of the instanton-sector path integral given in Appendix B, but now on the lattice. 
While to our knowledge this was never done in any detail, we stress that there is no 
conceptual difficulty here. In fact, the continuum change-of-variables (Appendix B.2) 
is a formal manipulation, because so is the original path-integral measure in terms of 
position eigenstates. In contrast, the corresponding change-of-variables on the lattice 
should be completely well defined (albeit technically more involved). 
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We recall that, in any non-perturbative sector, the collective coordinates are de- 
termined by listing the (exact or approximate) bosonic zero modes, that must be re- 
moved from the fluctuations spectrum to avoid infra-red divergences. The instanton- 
sector path integral on the lattice will, obviously, feature the same set of collective 
coordinates as in the continuum treatment. The anomalous SUSY Ward identity, 
eq. ( |3.10 ), will thus arise also within the above-sketched analytic lattice treatment of 



the instanton sector. (Of course, this applies also to the more general result, eq. ( |1.3| ).) 

It is interesting to re-examine the ambiguity of the previous subsection in a lattice 
context. On the lattice the instanton action is 

Statt (P) ~ ~YT~T\ + -2 Ma) , ( 4 -24) 
9 \P l 9o 

where go is the bare lattice coupling. The function f(p/a) accounts for the lattice- 
induced change in the instanton action when its size p becomes very small. First, 
identifying e of the previous subsections with the lattice cutoff a, the above advo- 
cated e/p — > prescription means that the lattice action should be chosen such that 
f(p/a) > 1 for p/a ~ 1 [|17[]. This condition implies that the minimal instanton size, 
while being vanishingly small in physical units, is infinitely large in lattice units in 
the continuum limit. (Typically, f(p/a) is a polynomial in a/ p. If /(p/a) ~ (a/p) 2n , 
unsuppressed instantons will be ones with p > agQ l ^ n .) 

The opposite prescription, namely e/p — > oo, corresponds on the lattice to choos- 
ing f(p/a) to be negative for p/a ~ 1. This means that the probability of finding 
lattice-size instantons is enhanced. In numerical simulations this enhances lattice arte- 
facts in physical quantities, which is undesirable. More seriously, this may become a 
problem of principles as we now explain. 

As discussed earlier, the matrix elements of the point-split operator Ojf gener- 
ically contain a piece that behaves like a delta-function 5(p). Suppose now that 
lattice-size instantons are not suppressed. When their Boltzmann weight is multi- 
plied by the effective S(p), a non-zero contribution to the Ward identity involving 
Ojf is obtained even in the (would-be) continuum limit g — > 0. As before, this 
contribution depends on the precise definition of Ojf but, moreover, it now violates 
Lorentz invariance. The reason is that it originates from lattice-size instantons which 
must be sensitive to the orientation of the lattice axes. (As a by-product, the SUSY 
Ward identity for Ojf will not be recovered on the lattice either.) 

In fact, the Boltzmann weight of lattice-size instantons must either diverge or 
vanish in the limit go — > 0. (It cannot stay finite without "infinite fine-tuning.") We 
conclude that, as a matter of principle, when the ambiguity of subsection 4.3 arises, 
there will exist a consistent continuum limit only if the lattice action is chosen such 
that lattice-size instantons are suppressed. As in the previous continuum treatment, 
the Ward identities for % f and Oz will then agree in the continuum limit of the 
lattice theory, and will both be anomalous. 
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5. Conclusions 



In this paper we have derived a general expression for the anomalous term in 
SUSY Ward identities (eq. (|1.3[)). We have analyzed in detail Ward identities in the 
one-instanton sector of SU(2) SUSY theories. We have found non-zero anomalous 
terms, which moreover cannot be modified by subtractions, both in SYM (eq. (|3.10|) ) 
and in SUSY theories with matter (eq. ( |3.32|) ). 

The SUSY anomaly arises from a Hilbert-space surface term, unlike the chiral 
anomaly which arises from non-invariance of the path-integral measure. In the SUSY 
case there seem to be no analog of the familiar anomaly-cancelation mechanism of 
the chiral anomaly. Similar anomalous Ward identities should exist with an SU(N) 
gauge group (subsection 3.4), and with any matter content. This would imply that 
the SUSY anomaly occurs in every asymptotically-free four- dimensional theory. 

Within a fully non-perturbative regularization (such as a lattice cutoff) the op- 
erator d^Sfn is necessarily non-zero. In the continuum limit, all matrix elements of 
d^Sfj, vanish to all orders in perturbation theory but, according to our results, d^S^ 
has non-perturbative matrix elements which are not zero. As a crucial check, this 
should be confirmed by a calculation which starts directly from the non-zero d^S^ of 
the regularized theory. 

Locality of d^S^ in the regularized theory implies its locality in the continuum 
limit. The local version of the anomalous Ward identity (eq. (|3.38| )) is consistent 
with this requirement. In eq. ( |3.3g| ) the correlation between the points and z must 
be mediated by a bosonic state. The ^-dependence, namely z u /(z 2 ) 2 , is that of (the 
derivative of) a massless boson propagator. However, the apparent single-particle 
propagation could also be due to two colinear particles. Usually the phase space for 
colinear propagation is zero, but in the triangle diagram the fermion and the anti- 



fermion do in fact become colinear in a special kinematic limit 22 . Whether a similar 



phenomenon takes place in the present case is another important question. Finally, 
the implications of this new anomaly on the physical spectrum of SUSY theories 
should be studied. 
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A. Notation 



For bosons we write: 

Bi(x) = b I (x)+j3 I {x), (A.l) 

where b;(x) and (3; denote respectively the classical and quantum parts of each Bose 
field Bj(x). We use a real notation where the generic index / runs over the gauge 
field and over both scalar fields $ and their complex conjugates $t. One has 

Bj & An, $ , $ t , 

bj ^> dp, ip, ^, (A. 2) 

Pi <^> a M , 0, f , 

where $ stands for all scalar fields in the theory. For the fermions we employ a 
Majorana-like notation |1| which is valid in vector-like SUSY theories, such as SQCD. 
For each pair of Weyl fermions in complex conjugate representations ip L , ip^ we define 

■0 = (-?/> L , <0 C = (-0^ ^ L ). We then write 

* f (s) ^ A , V , V> c , (A.3) 

where A is the (Majorana field) gaugino. We will also use the notation 

Qi(x) = bj(x) + qj(x) (A.4) 

where qt(x) stands for any quantum field, boson or fermion, and with the under- 
standing that bi(x) = for fermions (since there is no classical fermion field). The 
notation if?(x) will be used as a shorthand for the mode expansion of the fermion 
fields, cf. eq. 

We introduce a real inner product, defined for given values B^ and B^ of the 
bosonic fields as 

(bP\bP) = Jd 4 x (AWAW+J2($ {1)t ® (2) +® (2)t ® {1) )) ■ (A.5) 
The sum is over all scalar fields. For fermions the inner product is 

= J d 4 x (\V\M + Y,{^ 1) ^ 2) +^ 2) 4 1) )) ■ ( A -6) 

Here A = A T C con j etc., where C con j is the charge-conjugation matrix. Finally, the 
inner product of two adjoint-representation fields (e.g. the last term in eq. Q2.7|) ) is 
simply the integral of their product. 

The SUSY variation of any boson is linear in the fermions. We write this as 
follows 

<f a £f = r a fj*j. (A.7) 

The index a runs over all supersymmetries (i.e. we do not distinguish here between 
5 and 5). 
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B. Transverse-field Feynman rules 



In this appendix we derive Feynman rules for any non-perturbative sector of 
a weakly-coupled, renormalizable gauge theory. We restrict the discussion to the 
background Landau gauge. Both gauge and scalar classical fields will be considered. 
We face two technical problems. First, we are dealing with a coupled quantization 
problem, involving both (discrete) collective coordinates and (infinitely many) gauge 
degrees of freedom. This was discussed in the literature before, see e.g. |1| |2"0fl . The 
second complication arises because we do not assume that the background field is 
a solution of the classical field equations. Many relations used in the quantization 
around an exact solution do not apply here. 



The below Feynman rules do not involve constraints |]2T|. As a result there are 



linear (classical tadpole) terms in the action. The latter can be treated perturbatively 
provided the background field is an approximate solution of the classical equations. 
(Verification of the last statement must be done on a case by case basis. We will 
discuss the one-instanton sector of SUSY-Higgs theories elsewhere.) 

This Appendix is organized as follows. In subsection B.l we define the background 
gauge. Quantization in a given non-perturbative sector is discussed in subsection B.2. 
The mode expansion of the quantum fields, and various commutator formulae, are 
given in subsection B.3. In subsection B.4 we give a compact formula for the path in- 
tegral, which treats the ordinary collective coordinates and the infinitely-many gauge 
degrees of freedom on a uniform basis. This formula provides the basis for the discus- 
sion in Sec. 2. In subsection B.5 we defined the ghost action, and in subsection B.6 
we complete the Feynman rules by giving explicit formulae for the functional de- 
terminants and the tree-level propagators. Finally, in subsection B.7 we explain a 
complication that arises when one attempts to go to a general £-gauge. 



B.l. Background gauge 

We work in the background gauge whose explicit form will be repeatedly used. 
We first record the action of a gauge transformation parametrized by u a (x) on the 
classical field(s) 

6 u <p A = -igT A a B ip B u\ 5^ A = ig<fiT B a A u; a , (B.l) 

where 

Dp = <9 M — iga°^T a , (B.2) 
is the background covariant derivative. In compact notation this reads 

<J w 6 f (x) = n?(6(x))w°(x), (B.3) 

where fi~ is a linear differential operator that maps the Lie-algebra valued function 
uj a into the space of all (classical) bosonic fields. For each target Bose field, the action 
of can be read off from eq. (|B.1|) . 
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The quantum fields transform homogeneously. In particular, a c (x) transforms in 
the adjoint representation, with Tg c = —if a bc- We summarize all the transformation 
rules by 

6 U Q = (n a - igT a q)u a . (B.4) 
The background gauge is defined by 

nt°& = Dfa\ + ig £(^^0 - ^T\) = . (B.5) 

The classical field depends on the collective coordinates. In general, the ordinary ( n - 
derivatives of the classical field & n = db/d( n do not obey the background gauge. We 
introduce covariant (^-derivatives by making a compensating gauge transformation 

b, n = b n + Q a u a n , (B.6) 

where is defined by imposing the background gauge condition(s) 

& a b. n = 0. (B.7) 

A solution uj^ always exists, since 

Lf h = & a Q b = - (D 2 ) ab + ^{T a , T b } V , (B.8) 

is a positive operator. Under a gauge transformation, b- n transforms homogeneously. 
We define second (or higher) covariant derivatives of the classical field via e.g. 

b,n; m = b, n , m -igT a b, n u a m . (B.9) 



B.2. Quantization 

Formally, the path integral is defined in terms of position eigenstates, and has the 
measure VB(x)V^/(x). Quantization in a non-perturbative sector means a change 
of variables at the price of introducing a jacobian. The new variables include the 
collective coordinates ( n , the amplitudes of the quantum modes, and the gauge degrees 
of freedom u a (x) which decouple later. The change of variables is facilitated by 
introducing the following identity into the path integral 

1 = J d n C J Vuj{x) J 

xH5[( u b., m \B^- u b)} (B.10) 

m 

xJI 5[ u rt a (B^\x) - u b{x)) . 

x,a 

Here B^\x) denotes a finite gauge transformation of the bosonic fields generated 
by exp(— igu a (x)). The role of the first delta- function is to remove the infinitesimal 
variations b- n from the fluctuations spectrum, and to trade them with collective coor- 
dinates. The second delta-function fixes the gauge, as in the standard Faddeev-Popov 
procedure. 
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A frequently encountered technical problem is that one obtains ordinary ( n - 
derivatives in the Jacobian, and their replacement by covariant (^-derivatives requires 



some justification |Tj|. We overcome this problem by considering = u^(x;Q as 
a connection on the space spanned by the collective coordinates. We introduce a 
parallel transporter on this space 



U(x;C) =Vexp i-ig f ds n u a n T a \ , (B.ll] 



where V denotes path ordering and the integration runs along some path from the 
origin to (. In eq. (|B.10|) u b(x) denotes the action of the finite gauge transformation 
defined by U on b{x), and u ^ a = Q) a ( u b(x)). 

In the computation of the jacobian J, one must remember that B(x) and u a (x) 
are independent of the collective coordinates. The ^-dependence enters through 
b = b(x;() and U = U(x;Q. Thanks to the presence of U we obtain covariant 
(^-derivatives everywhere. One has 

U n = -igUT a LU*(() + path dependent terms , (B.12) 

where the first term on the r.h.s. comes from the variation of the end point. By choos- 



ing a sequence of paths judicious, the path dependent terms vanish in the limit 
After computing the (^-derivatives we set — u b = u (3. All the IA factors, as 
well as the explicit ci; a (x)-dependence, drop out now. The generalized Faddeev-Popov 
jacobian is 

, n ,( -Cmn -zgbl(y)T b P(y) \ 

\igbl(x)T a (3(x) ft^(x)(ft b (x)-igT b f3(x))S 4 (x-y) ) 1 ' 



where C = C + C\ and 



(Co)mn = (b-,m\b;n), (B.14a) 

{Cl) mn = (b ]m \P;n) = -(b ]m ;n\P) ■ (B.14b) 



In the derivation of eq. ( |B.13| ) we have also used eq. ( |B.7| ) and (fi a ) ;m = —igT a b- m . 
The matrix C is symmetric. For Co this is obvious. For C\ this follows using the 
commutation relations of the next subsection and eq. flB.5|) . 

We may now evaluate the delta-functions in eq. (p,10| ). The first delta-function 
removes the b- n modes from the bosonic fluctuations, yielding a factor of det _5 (Co). 
Similarly, the second (functional) delta-function removes the longitudinal modes and 
yields a factor of Det~2 (L gh ). The partition function thus reads 

Z= [ d n C Vq 1 J — 1 exp(-S) , (B.15) 

J det2(C )Det5(L gh ) 

where T>q = T>f3'D^. The quantum bosonic field obeys the gauge-fixing condi- 
tion ( |B.5| ) and 

(h rn \(3)=0. (B.16) 
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B.3. Mode expansion and commutators 



We now turn to a more detailed discussion of the functional measure. T>f3T>^> 
is defined by expanding each boson and fermion quantum-field in eigenmodes of the 
corresponding small fluctuations operator in the classical background. Repeating 
eq. (|2.1|) one has 

where q p denotes the amplitude of any quantum mode, and Xp( x ) is the corresponding 
eigenfunction. In the case of fermions, the mode expansion includes both the zero 
modes and the continuum modes. In the case of bosons there are usually only contin- 
uum modes. (However, the formalism can accommodate discrete bosonic modes too, 
should any survive after the removal of the b- n (x) modes.) 

Anticipating the need for ghost fields r] a (x) and f] a (x) (subsection B.5 below) we 
also give their mode expansion 



£c;(x) Vp , v a (x)=f / c a p (x)r ]p , (B.18) 



where the c p {x) are the (real) eigenfunctions of L^. The longitudinal bosonic modes 
which are eliminated by the background gauge are j3$(x) oc VL a Cp{x). 

The small fluctuations operators depend on the collective coordinates through the 
classical field. Through the eigenfunctions, this implies dependence of the quantum 
fields too on the collective coordinates. The ^-derivative of a quantum field is 



^X P> n{x)q p . (B.19) 



-p 

Covariant (^-derivative are defined by 

q, n = i n -igT a qu a nl (B.20) 

and similarly for the ghost field. 

We will also need various commutators. First, the commutator of two covariant 
derivatives of the classical field is 

b;m;n ^;n;m ^ ^mn i (B.21) 

where 

Fmn — U m,n ~ U n.,m + dfabc U m • (B.22) 

Including the commutator for any quantum field this reads (compare eq. ( |B.4j )) 

Q;m;n ~ Q;n;m = (0° ~ itfFftKn , (B.23) 

in agreement with our interpretation of io^ as the components of a connection in 
C-space. 
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Next we consider the commutator of a covariant ^-derivative and a gauge trans- 
formation. The local parameter of the gauge transformation is identified with a ghost 
eigenstate, u a (x) = c^{x) (see the next subsection). Since c£(x) is a function of the 
collective coordinates, one has 

(<LQ) ; „ " UQ;n) = (fi a - igT a q)F- n , (B.24) 

where 

Kn = M c = uj% - gf abc u a u b n . (B.25) 

(The above may also be applied to the improper gauge transformation that gener- 
ates the isospin zero modes, see Sec. 3.) Finally, for the commutator of two gauge 
transformations we have the familiar result 

S^Q - S^Q = (Q a - tgT a q)F^ 2 , (B.26) 

with 

K^=9f ahc ^\. (B.27) 

B.4. Uniform treatment of discrete and gauge collective coordinates 

Let us expand the local parameter of the gauge transformations u a (x) in terms 
of (its natural basis of) ghosts eigenstates 




(B.28) 



Each (unnormalized) longitudinal mode Q a Cp(x) is the infinitesimal variation of the 
classical field with respect to the gauge collective coordinate u p . This observation pro- 
vides the basis for a uniform treatment of all parts of the jacobian ( |B.13[ ). Taking the 
inner product of the continuous-index terms with ghosts eigenstates (which amounts 
to a unitary change of basis), the jacobian is rewritten as 

T-T)et( C ™ -w(b m \T b p\c») \ 

I (<$\SV*(Sl>-i { fl*0)\<i) ) { } 

Notice that (c£) ;n fi at /3 = by eq. (|B~5D . Hence the lower-left entry in eq. (|B~29p 
agrees with eq. (B.13| ). (We return to this observation in subsection B.7 below). The 



jacobian now takes the compact form 

J = DetC, (B.30) 

C M n — (Co)mn + (Ci)mn , (B.31) 

(C )mn = (b; M \b. N ) , {Ci) M n = -(P;m;n\I3) ■ (B.32) 



The capital indices m,n stand for the discrete indices m, n that label the ordinary 
collective coordinates, as well as for the indices p, q of the ghosts eigenstates, cf. 
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eq. ( B.2S ) above. (E.g., for the classical field one has explicitly b- M = b- m for m = m 
and b ;M = VL a c^ for m = p.) Substituting eq. ( |B.29| ) into eq. ( |B.15| ) we obtain 

d n CV q_^L e S, (B33) 
Det^ C 

Eq. ( |B.33| ) provides the basis for the general discussion of Sec. 2. 

B.5. Ghost action 

We complete the definition of the Feynman rules by giving expressions for the 
functional determinants and for the tree-level propagators. In this subsection we dis- 
cuss the ghost action. In a diagrammatic expansion, the discrete- and the continuous- 
index parts of the jacobian (p.l3j ) must be treated separately. Given a general square 
matrix 



one has det A = det 5ndet D where D = B 22 — B 2 \ Bx\B\ 2 (if B~{\ exists). Applying 
this formula to eq. ( |B.13| ) and introducing the ghost and the anti-ghost fields r\{x) 
and fj(x), the jacobian flB.13j ) can be written as 

J = det C J VrjVfj exp(-S fp ) , (B.35) 

where the ghost action is 

Sfy = J (lypl^lx,!,)^!/), (B.36) 

F ab (x, y) = F«\x, y) + f*(x, y) . (B.37) 
The tree-level and the interaction terms are 

^\x,y) = Lli(x)5\x-y), (B.38) 
?&(x,y) = -ign a ^x)T b p(x)5 4 (x-y) 

+g 2 tf m (x)T a p(x) C~i b\ n {y)T b p{y) , (B.39) 

where is defined in eq. ( |B.8| ). The ghost propagator G g h is defined by 

L gh G gh (x,y) = 6\x-y). (B.40) 

Next, the covariant derivatives b. m are 0(1/ g). Hence, a systematic expansion for 
C^ 1 (cf. eq. (p.39| )) is obtained by writing 

C- 1 =C 1 -C l C 1 C l + --- . (B.41) 

This may be represented in terms of discrete-ghost propagator and vertices, see Fig. 1. 
The Feynman rules depicted in Fig. 1 also generate the expansion of det (C/Cq). The 
partition function is now 



r det C 

\ d n ( VqVr]Vf] j j exp(-S - S ip ) . (B.42) 

J det5(C )Det5(L gh ) 
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B.6. Propagators, functional determinants and interactions 



The last step is to expand the bosonic fields in the action around the classical 
background. One obtains terms containing any number of quantum fields between 
zero (the classical action) and four. The interaction lagrangians are 

£f = \ SW h$ A* , (B.43) 

and 

C%(b, 0) = S?\b) P-, + l - S§f(b) P-Mk + y a $%i PtPsPkPl • (B.44) 

As explained earlier, the tadpole (linear) terms involving S^\b) = 5C B /5bi arise 
because we are not assuming that the background field is an exact solution of the 
classical equations. 

The bilinear terms serve to define the tree-level propagators and the functional 
determinants for bosons and fermions. The semi-classical measure (including the 
ghosts contribution) is 

exp(-Wi) = det5(Co) Det^(L gh ) DeH(L£) Det^(L F ) . (B.45) 

The one-half power of the fermionic determinant is due to our Majorana-like no- 
tation (Appendix A). As usual, exact zero modes are excluded from the fermionic 
determinant. Notice that, unlike in a general {-gauge, we have here the ghosts deter- 
minant to a one-half power. This is compensated by the absence of longitudinal-mode 
contributions in the transverse bosonic determinant. 

Bosons. We now turn to a detailed discussion of the propagators and the functional 
determinants, starting with the bosons. We focus on the propagator for those (gauge 
and scalar) fields that have a classical value. (For bosonic fields that do not have a 
classical part, the below derivation reduces to the textbook case, and the propagator 
is the inverse of the differential operator that occurs in the bilinear action.) Let the 
quadratic bosonic action be 

S%\b,P) = ±(p\L B \p). (B.46) 

The second-order differential operator L B = L B (b) depends on the classical field. 
Our task is to construct a bosonic propagator G B which is the inverse of L B on the 
transverse subspace. The latter is the complement of the subspace spanned by the 
b- n and by the infinitely-many longitudinal modes. It is convenient to consider first a 
general {-gauge, and then take the limit £ — > 0. This means that we first project out 
the finite-dimensional space spanned by the b- n for any £, and in the end remove the 
longitudinal modes by sending £ — > 0. 

We begin by considering the differential operator 

L% = L B + ^n a ^ a , (B.47) 
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where the last term is the contribution of the gauge-fixing action for £ 7^ 0. When the 
background field is not a classical solution, L B has no (exact) zero modes. Therefore 
it has an inverse Gb such that 

LiG B (x,y) = 6\x-y). (B.48) 

We will now express Gb in terms of Gb and b- n . (The below construction is similar 



to that of Levine and Yaffe pO ], see also ref. flgl .) We consider the bosonic gaus- 
sian integration in the presence of an external source K(x) at fixed values of the 
collective coordinates ( n - We introduce a bosonic field f3'(x) which obeys the orthog- 
onality conditions ( |B.16| ) but not the background gauge ( |B.5| ), as well as a completely 
unconstrained field f3"(x). We let 



)mn b n 



(B.49) 



denote the normalized covariant derivatives of the classical field. Now 



Zi{K) 



Vf3 f exp 



--(p'\Li\(3') + (K\P') 
1 



VP" Y[6 exp -UP"\L B \P") + (K\P») 



(B.50) 



Vf3" / da n exp 



kp^L B \p") + (K + i Y / a n b n \p") 



where, on the last row, we have introduced an integration over auxiliary variables a n 
to enforce the constraint. It is straightforward to perform the gaussian integration, 
first over /3'(x), and then over a n . The result is 



Z^K) = Det-i (L| x ) exp \-{k\G b \k) 
The bosonic determinant is 

Det (L B X ) = Det _1 (G B ) = Det (4) det (A" 1 ) , 

where 



(B.51) 



(B.52) 



A mn = (bm\G B \bn). (B.53) 

The bosonic determinant in eq. flETiSD is the £ -> limit of eq. (|B~52|) . (Det (4) 
contains an infinite power of the gauge- fixing parameter £; this power, which is for- 
mally equal to the number of longitudinal modes, cancels against a matching infinite 
power upon dividing by the free vacuum-sector determinant.) 

The normalization factor C 2 in eq. ( |B.49|) cancel out in the expression for the 
bosonic propagator, which reads 



Gb — Gb — Gb A Gb , 

A = Yl \ b ;m) A rnn (&;n| , 
mn 
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(B.54) 
^B.55) 



(A- 1 ) mn = {b;m\G B \b;n) ■ (B.56) 

At this stage the bosonic propagator obeys the orthogonality relation G B \b-n) = 
and the completeness relation 

4 G B (x, y)-J d A z P B {x, z)L% G B (z, y) = 6\x - y) - P B (x, y) , (B.57) 

where 

Ps(x,y) = ^~b n (x)~bi(y) 

n 

= Y,b i n^)(C 1 ) m nbl(y), (B.58) 

mn 

is the projector on the space spanned by b ]n . 

The last step is to take the Landau-gauge limit £ — > 0. A very useful identity is 

G sh rtL B G B = GghH 1 " - - StfG B , (B.59) 

which relates the longitudinal part of G B to the ghost propagator. Eq. (B.59| ) is 

derived by multiplying eq. (P-48|) on the left by G g h^ and integrating by parts. (If 
the classical field is an exact solution, the longitudinal modes are zero modes of the 
bosonic small fluctuations operator and the l.h.s. of eq. (P-59| ) is zero; eq. ( |B.59|) holds 
in fact for G B too, as one can check using eq. ( |B.54| ) and Q^A = 0.) Eliminating the 
gauge-fixing term from eq. (|B.48|) using eq. (|B.59|) , we rewrite the former as 

L B G B (x, y) — J d 4 z pH(x,z)L B G B (z,y) = 5\x - y) - P\x,y) , (B.60) 

where _ 

P^x,y) = nG gh (x,y)B, (B.61) 

is the longitudinal projector. Eq. ( [B.60| ) holds for any £, and has a smooth ^ — > 
limit. Finally, substituting eq. ( B.60|) into eq. (|B.54|) we find 

L B G B = l-KG B -P\l-L B G B ) 

= 1-(P B + P^(1-L B G B ). (B.62) 

The last expression exhibits again that the covariant derivatives b- n and the longitu- 
dinal modes have a similar role. This result is used in the diagrammatic identities 
of Appendix C.2. We comment that eqs. QB.571) , ( B-60| ) and ( p.62|) all have the 
generic form (1 —P)LG = 1—P, which is the completeness relation for a constrained 
propagator obeying P G = G P = 0. 

Fermions. In Majorana-like notation (cf. eq. ( A.6 )) the bilinear fermion action is 

SP(b,^>) = ^\L F \^). (B.63) 
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The euclidean Dirac operator Lp can be chosen to be hermitian [|18 |. 

In a generic SUSY theory that contains both explicit mass terms for matter fields 
and a scalar VEV, the fermion spectrum contains no exact zero modes, and Lp has an 
inverse Gp such that L F Gp(x,y) = 5 A (x — y). Alternatively, if there are no explicit 
mass terms, Lp anti-commutes with the generator of some (typically anomalous) R- 
symmetry even in the presence of a Higgs field. In a sector with a single classical 
object (one unit of topological charge) the classical field is spherically symmetric (or 
can be chosen to be so, in case it is not an exact solution). This implies the existence 
of a conserved index for each value of the total angular momentum. In this case the 
fermion spectrum contains exact zero modes. One has 

L F G F (x, y) = 5\x -y)- P^\x, y) , (B.64) 

where 

Pr ct (x, y) = £ X f°(x) Xi 01 (V) , (B-65) 

i 

is the projector on the exact fermionic zero modes xf°- 

In general, the fermion propagator Gp in the instanton sector contains a contri- 
bution from approximate zero modes too. The latter correspond to those zero modes 
of the exact instanton solution which couple through the explicit mass terms and/or 
the Higgs field. For practical calculations, it may be convenient to separate out the 
approximate zero modes. To this end one needs fermionic propagator and determi- 
nant which are constrained to the complement of the subspace spanned by both exact 
and approximate zero modes. These are defined via the same algebraic construction 
used above for bosons. 

For the general diagrammatic proof of SUSY Ward identities (Appendix C) it 
is convenient to use the propagator Gp. The projector on the exact fermionic zero 
modes p^ xact arises at intermediate steps in the derivation. However, the sum of all 
terms involving p^ xact is zero due to anti-symmetry. 

Recursion relations. A key role in the renormalization of (SUSY) Ward identities 
is played by recursion relations between the non-perturbative and free propagators. 
Let H = Ho + V be a Schrodinger operator. Assume that both Ho and H have no 
zero modes, and let Go and G be the corresponding Green functions. One has the 
textbook relation 

G = (l-GV)Go, (B.66) 

which may be iterated a finite number of times. In our case, eq. ( |B.66| ) may be applied 
to the ghost propagator (eq. ( |B.40|) ) and to the unconstrained bosonic propagator Gb 
(eq. ( P38D ) for £ ^ 0. 

The Born series is obtained by iterating eq. (|B.66| ) infinitely many times. How- 
ever, the series do not converge if Ho and H have different zero-modes spectra. Thus, 
a generalization of eq. ( P-66|) is needed for the fermion propagator Gp which obeys 
eq. ( |B.64j ). We first write 

L F = L™ c + V F , (B.67) 
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where L V F C stands for the free (vacuum-sector) Dirac operator(s) for all fields partic- 
ipating in the exact zero modes. Its inverse, the free fermion propagator, is denoted 



G v p c . We now claim that 



G F =(1- P£ xact - G F V F )G^ . (B.68) 



Let A denote the difference between the r.h.s. and the l.h.s. of the above equation. 
One has A(L F — V F ) = AL V F C = 0. This implies A = since L™ c has no zero 



modes. One can also verify directly that the r.h.s. of eq. ( |B.68| ) is a right-inverse: 



L F A = P| xact (l + V F Gf) = P F xact L F Gf c = 0. The role of the above recursion 
relations in the renormalization procedure is discussed in subsection 3.3. 

B.7. The case £ / 

As can be seen from the previous discussion, it is natural to work with a transverse 
bosonic field. Since eq. ( |1.3| ) deals with gauge-invariant operators, it should hold for 
any £, and one may want to check this explicitly. However, there is a technical 
difficulty that has prevented us from generalizing the Feynman rules to £ ^ 0. 

As explained below eq. (p.29| ), one must use the gauge condition (B.5|) in order 
to prove that eqs. ( |B.13| ) and ( |B.29| ) define the same jacobian. For £ ^ 0, however, 
the constraint £V(3 = is not respected. Therefore, eqs. (p.l3j ) and ( p.29|) do not 



define the same jacobian. We did not attempt to resolve the question of what is the 
correct jacobian for £ ^ 0. (The gauge-fixing constraint is used also in the calculation 
of 5u a , cf. eq. ( |2.8|) ; therefore it is likely that 5u a too will pick up terms proportional 
to for £ 7^ 0. In the calculation of eq. (|3.10| ) it is in fact possible to work with 



any £-gauge, since the above subtlety is relevant only at the next order.) 



C. General diagrammatic proof of SUSY Ward identities 

In this appendix we give a general diagrammatic proof of eq. (|1.3j ) . The underlying 
algebraic structure is similar to Sec. 2. However, the bilinear part of the action plays 
a special role, being the starting point for any diagrammatic expansion. As a result 
things are technically more involved. 

In more detail, besides algebraic manipulations, what entered the results of Sec. 2 
is the physical boundary conditions of the quantum fields. In a path-integral context 
we had to invoke a finite-volume cutoff in order to impose these boundary conditions. 
However, the propagators automatically obey correct boundary conditions in infinite 
volume. Therefore a diagrammatic expansion should be consistent with the vanishing 
of SS and Sfi without having to resort to a finite- volume cutoff. 

C.l. No gauge fields 

In this subsection we prove eq. (|1.3| ) in the physically less interesting case of 
(a non-perturbative sector of) a SUSY theory with no gauge fields. This simplifies 
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matters in a number of ways. First, there is no need to introduce a ghost field. In 
addition we now deal with ordinary (in place of covariant) derivatives with respect to 



the collective coordinates, which commute. Expectation values are given by eq. (1.2) . 



which results after performing the integration over the quantum fields in eq. ( |B.33j ) 



with C — > C (cf. eqs. (B.14) and ( |B.31|) ) . The semi-classical measure (eq. QB.45|) 
simplifies to 

exp(-Wi) = det^ (C ) DeH(L£) Det5(L F ) . (C.l) 

The interaction terms arising from the action are given by eqs. flB.43j ) and ( |B.44j) . 
Additional vertices arise from the expansion of det (C/Cq), cf. eq. ( |B.41| ). 



The basic strategy is to construct separately Ward identities for the various terms 
in the transformation laws ( ]2.1UQ and ( |2.11| ) of the quantum amplitudes. The final 



Ward identity is obtained by piecing together the individual Ward identities after 
suitable arithmetical manipulations. We begin by considering all diagrams that define 
the expectation value of some (multi)local operator O, times an insertion of 

(6Q\L\q) = (T$\L B \P) + (5V\L F \%) , (C.2) 

In each diagram we act with Lb {L f ) on the propagator whose one leg corresponds to 
(3 (¥) of eq. ( |C.2j ). Using eqs. ( |B.54| ) and ( |B.64|) we obtain the diagrammatic identity 

(5Q\L\q)0 = _p F exact ^2 -S f O^ + (r£|l - AG B \x) . (C.3) 

Here 

x = ?° + L!*c->-s B )o, (c.4, 

8(3 V d(3 J 

and S B and S F are defined in eq. ( |3.20|) . The = sign denotes equality under the expec- 
tation value defined by integrating over the quantum fields only. (This is an equality 
between disconnected diagrams, which correspond to exp(— W 2 ) (■ ■ in eq. ( |1.2j ).) 

As in Sec. 3, all terms with p^ xact in eq. (|C.3|) cancel each other by antisymmetry. 
Using SUSY invariance of the action, as expressed in eq. ( |3.19| ), we find 



SO 



(tv\aG b \x) - (vty^Otr ^C- 1 ^ . (C.5) 



Eq. ( p.5|) says that the fixed-C expectation value of 50 does not vanish because of 
terms arising from the jacobian in eq. (|B.33[), or from the constraints obeyed by the 



bosonic field (eq. ( |B.16| )), or from both. Our task will be to simplify the r.h.s. of 
eq. O. 

The hint comes from eqs. (|2.10|) and (|2.11 ). These equations show that, besides 



the expected piece arising from 5Q, the variation of a quantum amplitude contains 
additional pieces needed to compensate for the ^-dependence of the field. We thus 
proceed by writing down a diagrammatic identity for the variation 5\q(x) = 5( n q, n (x) 
(cf. eq. ( |B.19| )). To this end, we consider the new insertion 

^Kn{'{in\L\q)'. + :(q\L\i n ):} , (C.6) 
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where 5( n is given by eq. (2.9) with C — > C. The normal ordering prescription means 
that the q and q >n must not be contracted with each other. We momentarily introduce 
a finite volume cutoff, which allows us to drop total derivatives of currents involving 
g n as in Sec. 2. As before, we let L act on the propagator attached to the q-leg of the 
insertion, and the resulting terms with p^ xact cancel out. We obtain the diagrammatic 
identity 

= 6( n ^:(q\L n \q):-:SW:}o 

+*»(MI-^) + (M1P) < C7 » 

+ 5( n (/3 n \l-AG B \x) + (pJl-AG B \ ) O . 

VI I dp J 

In a finite volume one has 

where the (^-derivatives act on the (discrete) eigenvalues (compare eq. ( j2.20| )). Also 
recall that each propagator has the mode expansion 

G(x,y) = J2x P (x)\; 1 xl(y)- (C9) 

v 

This allows us to rewrite eq. ( |C.8|) as 
= l - :(q\L n \q): 5( n O 

d(6( n 



- (An I AG B \5( n X + (CIO) 



dp 

+ eMS' mt ) : ( |: | G, n | ^ ) : exp(-S int ) 6( n O 

Eq. ( J(J.1U| ) holds in the infinite- volume limit. On the third row, S mt contains all the 
interaction terms (that arise from the expansion of both the action and the jacobian). 
The functional operator on that row differentiates (the product of) all propagators in 
a given diagram with respect to ( n . On the last row, xf° denote the exact fermionic 
zero modes (see eq. ( B.65| )) and q[° are the corresponding Grassmann amplitudes. 



The role of the last row is to differentiate all (exact) fermionic zero modes with respect 

tO Cn- 

It may be helpful to examine the content of identity (|C.10| ) in a simpler case, 



where it is applied to a scalar field that does not have a classical part. The second 
and forth rows are then missing, and the identity reduces to the familiar relation 



G, n (x,y) + J d 4 zG(x,z)L, n G(z,y) = 0, (C.ll) 
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which holds in any volume. As for the terms on the second and forth rows of 
eq. ( |(J.10[ ), all of them involve localized sources which guarantee a smooth infinite 



volume limit. (At coinciding points, x = y, eq. ( U.ll|) may however be violated by 



the ultra-violet cutoff (see subsection 3.3). If angular momentum cutoff is employed, 
the partial-wave propagators obey eq. (|C.11J) even for x = y. In that case the l.h.s. 
of eq. ( |C.11| ) is formally <5 n (0), and the spectral trace in eq. Q2.33|) may be regarded 



as its integral.) 

Consider now the sum of eqs. (|C.5|) and (|C.10| ). In order to arrive at eq. ( |1.3|) we 



need, besides (-derivatives of the propagators, also (-derivatives of all elements of the 
Feynman rules that depend explicitly (not through the eigenmodes) on the classical 
field. The latter must come from the r.h.s. of eq. ( |U.5| ) plus the first two rows on the 
r.h.s. of eq. flC.lOQ . Verifying that this is indeed the case is essentially an arithmetical 
task. The reader who is not interested in these algebraic details is advised to skip 
to subsection C.2, where we explain how we deal with the extra complications of a 
gauge theory. 

We start with the terms involving X in eqs. ( |C.5| ) and ( |C.10| ). Writing AG B = 
(AG B - P B ) + P B (cf. eq. ( jBjp ) we get 

(r£ - 5( n p, n \AG B \x) = 5Cn (b n + d n \x) , (C.12) 

where 

d m = ((Co) m iAi n — 5 mn ) b >n (C.13) 

and Ai n is defined in eq. (|B.56| ). We now consider the r.h.s. of eq. (|C.12| ) term by 
term using eqs. (|R43|) , ( ggg) , (|3~20| ) and (£J). First, b n dO/dp = b >n dO/db 



since O = 0(b + [3,^). This gives the (-derivative of the classical part of O. 
Next {b hn \sf\b)) = b, n dSW/db gives the (-derivative of the classical action, and 
{bi,n\sfj kL ^j^Kh) gives b tn d/db acting on Sf~?' (b) Pi fliPk- 

Next consider those terms in S B arising from eq. (|B.43| ) and from the remain- 
ing (the middle) term in eq. ( p.44j) . Combined with the first term on the r.h.s. of 
eq. ( |CTLJ ) we find 

:(q\L n \q):- (bj, n \si%(b)qrq & ) = -Tr(-) F L n G. (C.14) 

We used Lja, n = 6.; >n 5 7 -j)-.(6). Including a one-half factor, this gives the logarithmic 
(■-derivative of the functional determinants, except for those bosonic fields that have 
a classical part, see below. 

We next turn to the d n term in eq. (|C.12j ) . In the limit where the classical field is 
an exact solution, d n vanishes. (This can be worked out using Appendix B.6.) In our 
Feynman rules, not having an exact classical solution is reflected in the tadpole terms 
and in the peculiar formulae for the bosonic propagator and determinant. We will 
see that the <i n -dependent terms are indeed related to those elements of the Feynman 
rules. First, similarly to eq. (|C.12| ) we rewrite the second term on the second row of 
eq. (|C.10| ) as O times 

-(P !n \AG B \d(5( n )/dp) = -(C^UiC^bmldiSC^/dp) 

+ (Smn ~ (C^Udn) (d m \d(5( n ) / d@) .(C.15) 
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The <i n -dependent term in eq. ( |C.12| ) plus the first term on the second row of eq. (|C.15| ) 
give 

(d n \K n X+ d -^do) = -(b >n L B G B \8C n X + ^-o) 



± -(b n dSP/db\fc)5( n O. (C.16) 



This is the ^-derivative of the classical source sP{b). The remaining <i n -dependent 
term in eq. ( |C.15| ) plus the r.h.s. of eq. ( |C14j ) (for those fields that have a classical 
part) give 

- - 5( n Tr L n G - {C, l ) ml C ln (d m \ d(5( n )/dp) = 5( n — log DetT* (L B ) . (C. 17) 

(Note that the definition of the bosonic determinant (|B.52|) involves the normalized 
modes of eq. ( |B^|) .) 

We are now left with three terms: the last term on the r.h.s. of eq. ( |C.5| ), O times 
the first term on the r.h.s. of eq. (|C.15|) , and 

6( n O (b, n | d -^h^C u } = -5( n O (b n \b Al ) Cg 1 , (C.18) 

which is the last term that has remained from the r.h.s. of eq. ( |C12 ). In order to 
complete the derivation of eq. (|1.3|) we must show that the sum of these three terms 
is equal to O times 

(b n \d(5( n )/db) +5CnCf k 1 (b n \dC kl /db) -^5( n (C l ) lk (b n \d(C ) kl /db) . (C.19) 

This expression amounts to b^djdb acting on (the logarithm of) the jacobian and 
on S( n itself. The needed result follows after a few algebraic manipulations using the 
commutativity of ordinary (^-derivatives. This completes the diagrammatic proof of 
eq. ( |1.3| ) for a SUSY theory that contains only scalar and fermion fields, and no gauge 
fields. 



C.2. Gauge theories 

Before we turn to the full-fledged diagrammatic proof of eq. (|f .3| ) in a gauge 
theory, there is one more exercise that we can do. In the discussion of the previ- 
ous subsection, let us make the (formal) replacement C — > C, m, n — > m,n etc (cf. 
Appendix B.4). In other words, we consider both the discrete and the continuous 
collective coordinates. Accordingly, we relax the assumption that Cjv-derivatives com- 
mute. Repeating the entire derivation we now arrive at eq. ( |1.3|) , but with one extra 
term which is left after the algebraic manipulations. This term is (O times) the 
commutator from eq. ( [2.33 ). 

The careful diagrammatic treatment of the boson and the fermion fields has con- 
firmed the vanishing of the corresponding spectral traces in eq. ( |2.33|) . This is simply 

46 



because no such terms are left in the final diagrammatic answer of subsection C.l. 
Now, as discussed in Sec. 2, the commutator is really a spectral trace over the ghost 
field. As long as we treat the continuous-index part of the jacobian formally, we 
cannot escape the appearance of this commutator. The full gauge-theory Feynman 
rules are needed in order to confirm the vanishing of the ghost-field spectral trace. 

The general strategy will be the same as in Appendix C.l. We begin by establish- 
ing a number of key diagrammatic identities which correspond to various pieces in the 
transformation rules of the quantum amplitudes (|2.10|) and (|2.11|) . The rest amounts 



to very lengthy arithmetical manipulations that we omit. The gauge-theory partition 
function is defined by eq. ( |B.42|) . The semi-classical measure is given by eq. ( |B.45|) . 



The interaction lagrangians arising from the action are eqs. ( |B.43| ) and ( |B.44j ). The 



ghosts lagrangian as well as the vertices arising from the discrete part of the jacobian 
are given in Appendix B.5. 

The first diagrammatic identity that we need is a generalization of eq. ( |C.5|) . 
Again we consider insertion (|C.2|) . Performing the same steps as in Appendix C.l 
and using eq. ( |B.62| ) we obtain 



<)0 = f L * I (F'i; Pi )( I L„ On) II X (?) 




cr 1 --^e>) , (c.20) 

where X is defined in eq. ( |U.4j ) and the ghost interaction action can be read 
off from eqs. ( |B.36|) and ( B.39|) . On top of the r.h.s. of eq. (|C.5D , we see here the 
contributions of the ghost action. 

The second diagrammatic identity corresponds to the variation Sl° v q(x) = S( n q ;n {x) 
Following the same steps as in Appendix C.l we obtain 

^ l - :(q\L, n \q):5CnO (C.21) 



- (An I (Pb + P«)(l - L B G B ) I SC n [x- d M-0)+ ^ O) 
+ exp(S int ) : ( ^ | G. n | ^) : exp(-5 int ) S( n O 



As with eq. ( |C.1(J| ), the functional operator on the third row differentiates all bosonic 
and fermionic propagators with respect to ( n , while the last row does this for the 
exact fermionic zero modes. For the (^-derivative of the ghost propagator we may use 
(the covariant version of) the simpler relation eq. ( |C.11| ). 

The last two identities correspond to the variation 5 2 q(x) = —ig5u a (x)T a q(x). 
This variation looks like a gauge transformation of the quantum fields, except that we 
have replaced the local parameter u a (x) by its SUSY variation 5u a (x), see eq. (|2.8|) . 
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For bosons the identity reads 



= y :(5u a T a p\L B \p): O 



( 



- ig ( 5uo a T a (3 



1-{P B +PW){1-L B G 



B 



X 



as 



mt 
fp 



dp 



O 



) 



(C.22) 



+ | d 4 * (P B (x,2)+pll(x,2)) L B G B (z,y) 



S 4 (x — y) — Pb(x,v) - P ll (x,y) 



O. 



and for fermions 



= j :(5uj a T a y\L F \y): O - tg (^5uj a n ^-S f O^ 



, 4 d(6u a (x)) 
iq / a x 



T a ^(x) O. 



(C.23) 



With the above diagrammatic identities at hand, eq. ( |1.3| ) follows after lengthy 
arithmetical manipulations which will not be repeated here. Apart from similar steps 
to those of Appendix C.l, the commutator formulae (Appendix B.3) play an impor- 
tant role. In addition we use gauge invariance of the theory in general, and of the 
operator O and the functional determinants in particular. (The detailed discussion 
of Sec. 3.2 provides a leading-order explanation on how eq. (|1.3|) works in a gauge 
theory.) 



D. Translation invariance 



It is instructive to see how translation invariance is proved using the general 
formalism of Sec. 2. Let us define S^Q = Q ;A1 . Using eq. ( |2.2j ) one easily finds that the 
translation collective coordinates Xq transform according to S^x^ = —5^. Moreover, 
all other collective and quantum variables are invariant. The corresponding Ward 
identity (cf. eq. (|2.15|) ) for a gauge-invariant local operator reads 



d'xo^iO) 



8xq 



' x 







(D.l) 



We see that the Hilbert-space surface term coincides with a spacetime surface term, 
which vanishes thanks to locality (and gauge invariance) of O. 

Let us now examine the local Ward identity involving the conserved energy- 
momentum tensor T^ u (compare subsection 3.5). Define S Zfl Q(x) = 5 A (x — z)Q ;fl (x). 
Now other collective coordinates besides Xq transform as well. In particular, in the 
one-instanton sector the scale collective coordinate p transforms to leading order as 



g 2 da c v {z) 



167r 2 dp 
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(D.2) 



Using the explicit expressions for the instanton field it easily follows that S^p(z) is a 
total ,2-derivative, whose integral is zero. Suppose now that, for a given operator O, 
the p — > surface term involving 8^p{z) is finite. Being a total spacetime derivative, 
this term can be absorbed into a redefinition of the matrix element of T^ u (z). (An 
additional total spacetime derivative which is absorbed into the matrix element of 
Tfiu(z) comes from the local variation of the measure.) 
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H — I — I — I — h 



fermion (gaugino) propagator 
boson (gauge field) propagator 
ghost propagator 

any (boson, fermion or ghost) propagator 

C o (or "discrete ghost" propagator, cf. eq. (B.14a)) 



G 

D- 



fermionic zero mode 



classical field 



(unnormalized) bosonic zero mode 

the differential operator Q. (cf. eq. (B.3)) 



last term on the r.h.s. of eq. (3.21) 
eq. (3.25) 



v 



d (eq.(B.14b)) 




off-diagonal term of the jacobian (B.13) 



second term on the r.h.s. of eq. (3.26) 



Figure 1: Feynman rules. 



51 



Figure 2: Trivially vanishing leading order. The three fermionic zero modes coming 
together saturate the operator Ojf. The other two external legs correspond (on the 
l.h.s.) to the operator AA or (on the r.h.s.) to its leading-order SUSY variation. In 
all figures, the indicated equalities hold after antisymmetrization (see fig. 4). 







Figure 3: Bubble diagrams. The first and last diagrams contain a sum over 
boson, fermion and ghost contributions. The third diagram involves the off-diagonal 
terms of the jacobian eq. ( |B.13j ). The last diagram involves (one-loop) self-energy 
counterterms. 
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Figure 4: Antisymmetrization. The fermionic projectors (eq. ( |B.64j ) or eq. ( |3.7j )) 
cancel each other after antisymmetrizing with respect to the zero modes and the 
insertion of d^Sj®\ 





Figure 5: A diagram where d^S^ is connected to a fermionic interaction vertex. 
Here and in the following figures, those external legs which are simply saturated by 
fermionic zero modes are not shown. As before, the equalities hold after antisym- 
metrization (including those zero modes which are not shown here, see Fig. 4). The 
filled circle (here in the middle diagram) indicates that a bosonic field is attached to a 
fermionic leg of a vertex. The index-matching differential operators (or matrices) can 
be read off from the appropriate terms in eq. ( |3.22| ). On the r.h.s. this is traded with 
an insertion of eq. ( |3.25| ). These diagrams contain a loop (the almost-closed triangle) 
provided the two external legs are at the same spacetime point. 
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Figure 6: The third diagram on the r.h.s. is a first contribution to the p-derivative 
of ((6 a p) Ojf(0) XX(y)) p . In that diagram, the two rectangles connected by the think 
dashed line constitute the bosonic projector Pb, cf. eqs. ( ^78|) and ( p3.58|) . The first 
diagram on the r.h.s. is a contribution to 50, while the second diagram vanishes by 
a Fiertz rearrangement (if the need for regularization is ignored). The last diagram 
involves the longitudinal projector ( |B.61[ ), and together with other diagrams add up 
to zero by gauge invariance, see text. 






r0\ 



--{}*- 




Figure 7: Another diagram with an insertion of eq. ( |3.25| ). Again, the second 
diagram on the r.h.s. vanishes after a Fiertz rearrangement, the third diagram is a 
contribution to the p-derivative of ((S a p) Ojf(0) XX(y)) p , and the last diagram cancels 
against other diagrams by gauge invariance. 
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Figure 8: One-loop tadpole diagrams. The loop stands for the sum of boson, fermion 
and ghost contributions. The first diagram on the r.h.s. is a contribution to SO. In 
general, the third diagram gives the (^-derivative of the logarithm of the functional 
determinants. In SYM the functional determinants cancel each other exactly, and 
accordingly the one-loop tadpole is zero. In the second diagram on the r.h.s., it is 
convenient for our purpose to consider the contribution of each field separately. The 
fermion loop is by itself zero after a Fiertz rearrangement. The boson loop appears 
in Fig. 9. The ghost loop cancels against a sum of diagrams with the longitudinal 
projector, see text. 




Figure 9: The diagram on the r.h.s. is (which makes a contribution to 50). 

The first diagram on the l.h.s. is from Fig. 7, and the second one from Fig. 8. 
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Figure 10: The third diagram gives the p-derivative of (the classical field contained 
in) a covariant derivative D^. The first two diagrams contribute to 50. 




Figure 11: Diagrams containing the matrix C% (eq. ( |B.14b| )) from the discrete- index 
part of the jacobian. The first diagram on the r.h.s. is a contribution to SO. The 
second diagram gives the p-derivative of the bosonic zero mode contained in Sp itself. 
The third diagram gives the p-derivative of logdet^Co (see text). 
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(b) x o X 



Figure 12: (a) A fermion self-energy diagram (cf. Fig. 7). (b) The corresponding 
counterterm diagram. 




(a) (b) 



Figure 13: (a) A diagram containing a counterterm for the operator AA. (b) A 
counterterm for the operator 5 (XX) = (i/y/2) o^F^ X. The corresponding one-loop 
diagrams are in Figs. 4, 5 and 6. 
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